MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


INVESTIGATION  OF  PHYSICS  OF  SYNTHETIC 
APERTURE  RADAR  IN  OCEAN  REMOTE  SENSING 

TOWARD  84/86  FIELD  EXPERIMENT 


♦ 


INTERIM  REPORT 
VOLUME  II: 

CONTRIBUTIONS 
OF  INDIVIDUAL 
INVESTIGATORS 

MAY  1986 


Approved  for  Public  Release;  Distribution  Unlimited 


"TOWARD  84/86"*  FIELD  EXPERIMENT 
INVESTIGATION  OF  PHYSICS  OF  SYNTHETIC  APERTURE  RADAR  IN 
OCEAN  REMOTE  SENSING 


INTERIM  REPORT 

VOLUME  II:  CONTRIBUTIONS  OF  INDIVIDUAL  INVESTIGATORS 


OMAR  H.  SHEMDIN, 
PRINCIPAL  INVESTIGATOR 


LIST  OF  PARTICIPANTS: 


ERIM: 

D. 

Lyzenga 

NRL: 

W.  Garrett 
W.  Keller 

JPL: 

W. 

Brown 

W.  Plant 

D. 

Held 

D. 

Hoff 

ORE: 

P.  Hwang 

T. 

W.  Thompson 

E.  Tajirian 
D.  Thomas 

NPGS: 

K. 

Davidson 

G. 

Geernaert 

SIO: 

R.  Guza 

NOSC: 

P. 

Hanson 

UKa: 

H.  Chaudry 

J. 

Wright 

R.  Moore 

UMd: 

R.  Harger 

FOR 


OCEAN  SCIENCES  DIVISION 
OFFICE  OF  THE  CHIEF  OF  NAVAL  RESEARCH 
ARLINGTON,  VIRGINIA 

(Technical  Monitor:  Hans  Dolezalek) 


MAY  1986 


(TOWER  OCEAN  WAVE  AND  RADAR  DEPENDENCE) 


TABLE  OF  CONTENTS  (  ^ 

VOLUME  I:  Pag 

EXECUTIVE  SUMMARY 

I.  INTRODUCTION  .  1 

II.  EXPERIMENTAL  STRATEGY  AND  MEASUREMENTS  .  4 

III.  DATA  ANALYSIS  APPROACH  .  23 

IV.  IN-SITU  MEASUREMENTS:  OBJECTIVES,  DATA  SUMMARY 

AND  EARLY  RESULTS  .  28 

V.  RADAR  BACKSCATTER  MEASUREMENTS:  OBJECTIVES, 

DATA  SUMMARY  AND  EARLY  RESULTS  .  58 

VI.  SAR  IMAGING  ARGUMENTS . 73 

VII.  SAR  MEASUREMENTS:  DATA  SUMMARY  AND  EARLY  RESULTS  .  86 

VIII.  SUMMARY  AND  PRELIMINARY  CONCLUSIONS  .  105 

IX.  ANALYSIS  PLANS  AND  A  PERSPECTIVE  ON 

FUTURE  RESEARCH  .  107 

REFERENCES 


APPENDIX:  DATA  LOGS  AND  ENVIRONMENTAL  DATA  IN  SUPPORT  OF  MATERIAL  IN  TEXT 
VOLUME  II: 

INTERIM  REPORTS  OF  INDIVIDUAL  INVESTIGATORS 

I.  ^  SURFACE  GRAVITY  WAVE  MEASUREMENTS,  by  R.  T.  Guza 

II.  c  DEVELOPMENT  AND  UTILIZATION  OF  A  SURFACE  ENERGY  MEASUREMENT 

SYSTEM  IN  TOWARD  84/85^ by  W.  D.  Garrett 

III.  WAVE  FOLLOWER  MEASUREMENTS  DURING  TOWARD  84/85,  by  P.  A.  Hwang 

IV.  f  THE  PROPAGATION  OF  SHORT  SURFACE  WAVES  ON  LONGER  GRAVITY  WAVES, 

by  M.  S.  Longuett-Higgins 

V.  TOWARD  METEOROLOGY  MEASUREMENTS ,  by  G.  L.  Geernaert 

VI.  PROCESSING  OF  JPL  SAR  FRAME:  AZIMUTHAL  WAVES  ON  31  OCTOBER  1984,. 

by  E.  Tajirian  -y 

VII.  ^SAR  IMAGERY  SIMULATED  FROM  TWO-SCALE  RADAR  WAVE  PROBE  RETURN,  . 

by  W.  E.  Plant 

VIII.  V SAMPLE  PREDICTIONS  AND  SIMULATIONS  OF  SAR  OCEAN  IMAGERY  *  / 

by  R.  0.  Harger 


APPENDIX  I 


Surface  Gravity  Wave  Measurements 


R  .  T  .  G  u  z  a 

Scripps  Instiution  of  Oceanography 
Center  for  Coastal  Studies 
Code  A009 

La  Jolla,  CA  92093 


The  characteristics  of  surface  gravity  waves  within  the 
California  Bight  are  profoundly  influenced  by  the  topographic 
setting.  The  Channel  Islands  and  the  orientation  of  the 
coastline  shelter  the  coast  from  many  wave  directions  (Figure  1). 
Pawka  [1982,  1983]  and  Pawka  et  al  [1984]  discuss  extensive 
measurements  of  the  directional  wave  climate  at  Torrey  Pines 
Beach,  approximately  10  km  north  of  the  NOSC  tower  (Figure  1). 

The  waves  at  these  two  sites  are  similar,  so  the  work  at  Torrey 
Pines  is  relevant  to  the  TOWARD  experiments.  At  Torrey  Pines, 
only  a  narrow  directional  sector  (282°-290°)  is  directly  open  to 
the  north  Pacific  swel.l  which  is  an  important  source  of  coastal 
wave  energy  (Figure  1).  However,  the  presence  of  shallow  shoals 
and  banks  in  the  island  vicinity  (notably  Cortez  and  Tanner  Banks 
in  Figure  2)  theoretically  provide  a  refractive  energy  source. 
Relatively  high  angle  north  swell  is  refracted  by  the  banks 
towards  Torrey  Pines  (Figure  2).  Thus,  unimodal,  deep  ocean 
north  swell  theoretically  results  in  a  bimodal  directional 
spectrum  at  Torrey  Pines;  the  open  window  at  about  290°  and  a 
refractive  peak  around  255°.  Extensive  measurements  by  Pawka  and 
collegues  at  Torrey  Pines  have  firmly  established  the  existence 
of  the  southern  refractive  peak  (Figure  3).  These  data  were 
collected  in  March,  1977  with  a  400  m  long  multi-element  linear 
pressure  sensor  array  and  analyzed  with  high  resolution  data 
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adaptive  estimators. 

Similar  high  resolution  measurements  have  not  been  attempted 
at  the  TOWARD  site.  However,  theoretical  refraction  calculations 
for  the  Mission  Bay  entrance  Channel  (i.e.  the  NOSC  tower) 
suggest  a  comparable  bimodal  response  to  north  swell  (Figure  4). 
Note  that  the  angular  location  of  coastal  directional  peaks  is 
somewhat  insensitive  to  the  deep  ocean  direction.  North  swell 
energy  can  only  reach  coastal  waters  through  the  open  window  or 
via  banks  refraction. 

The  NOSC  tower  site  is  not  locally  sheltered  from  southern 
swell.  Thus,  although  perhaps  less  energetic  than  northern  swell 
outside  the  Channel  Islands,  southern  swell  can  be  an  important 
component  inside  the  Channel  Islands.  The  characteristics  of 
southern  swell,  in  the  northern  hemisphere,  are  not  well  known. 

In  a  classic  series  of  papers,  Munk  and  co-workers  established 
the  source  as  high  latitude  storms  in  the  south  Pacific  and 
Indian  Oceans.  The  superposition  of  a  typical  southern  swell 
directional  peak  (230°)  with  the  bimodal  north  swell  response 
results  in  a  complex  trimodal  directional  spectrum  at  Mission  Bay 
(Figure  5).  These  peaks  are  separated  from  each  other  by  about 
30°  in  local  deep  water  (i.e.  a  few  km  offshore).  Refraction 
reduces  the  angular  separation  at  the  tower.  A  long  linear  array 
would  be  required  to  resolve  the  directional  peaks  (Figure  5). 

Cost  constraints,  and  the  relatively  low  priority  given  to 
directional  long  wave  measurements  in  this  experiment,  did  not 
permit  the  installation  of  a  long  linear  array  in  the  TOWARD 
experiment.  Four  pressure  sensors  were  installed  on  the  tower 


legs,  about  7  below  the  surface.  These  sensors  were  used  to 
estimate  two  orthogonal  components  of  sea  surface  slope,  similar 
to  a  pitch  and  roll  buoy.  This  type  of  directional  measurement 
has  relatively  low  resolution,  less  than  the  150  m  long  linear 
array  shown  in  Figure  4.  The  resolution  of  pitch  and  roll  type 
systems,  and  a  variety  of  analysis  techniques  are  discussed  by 
Hasslemann  et  al  [1980],  Lawson  and  Long  [1983]  and  Oltman-Shay 
and  Guza  [1984].  Figure  6  shows  that  bimodal  spectra,  with  peaks 
separated  by  90°,  can  be  resolved  by  the  estimators  discussed  by 
Oltman-Shay  and  Guza  [1984].  For  angular  separations  less  than 
about  60°  the  peaks  are  not  resolved  and  the  angle  estimated  to 
have  the  maximum  power  is  near  the  weighted  average  of  the 
unresolved  peaks  (Figure  7) .  The  estimator  of  Hasselmann  et  al 
[1980]  also  resolves  widely  spaced  peaks.  Neither  estimator  can 
resolve  three  peaks.  It  appears  that  present  methods  are 
incapable  of  resolving  the  trimodal  directional  spectra 
anticipated  at  the  experiment  site  (Figure  7).  The  Oltman-Shay 
and  Guza  [1984]  estimator  is  used  below. 

The  present  analysis  is  focused  on  four  overflight  days  of 
phase  1  (17  Oct,  31  Oct,  4  Nov,  7  Nov  1984)  and  8  overflight  days 
of  phase  2  (13,  18,  19,  20,  24,  25,  26,  27  March  1985).  The  time 
periods  analyzed,  significant  wave  heights  (obtained  from  depth 
corrected  pressure  sensors  about  7  m  below  the  surface),  mean 
currents  and  wind  speed  and  directions  (NOSC  sensors)  are  shown 
in  Tables  1  and  2.  Sea  Surface  elevation  spectra  are  shown  in 
Figure  8.  Most  of  the  spectra  are  peaked  around  .07  hz  (31  Oct 
is  an  exception  with  a  peak  around  .11  hz).  The  majority  of  the 
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spectra  are  rather  broad,  with  significant  energy  spanning  the 
.1-.3  hz  frequency  band.  Days  with  a  relatively  narrow  spectral 
peak  are  4,  7  November  and  18,  24  March. 

Directional  spectral  have  been  estimated  for  each  frequency 
band  in  Figure  8.  For  example.  Figure  9  shows  estimates  for  four 
frequency  bands  (width  .0039  hz)  on  31  October.  The  three  lowest 
bands  correspond  to  peaks  in  the  energy  spectrum  (vertical  arrows 
on  31  October  panel  in  Figure  8).  It  is  reemphasized  that  these 
directional  spectra  are  smoothed  (actually  convolved)  versions  of 
the  true  spectra.  Nevertheless,  the  directional  spreads  probably 
provide  at  least  a  relative  measure  of  the  spreading  in  each 
frequency  band. 

It  is  not  practical  to  present  the  hundreds  of  directional 
spectra  which  have  been  estimated.  The  angles  with  maximum 
power,  and  the  angular  spread  at  half  maximum  power  are  shown  in 
Figure  10.  In  general  the  low  frequencies  show  the  smallest 
directional  spread  reflecting  both  refractive  narrowing  (which  is 
stronger  for  long  waves)  and  the  well  known  tendency  for  long 
waves  to  have  narrower  beam  even  in  deep  water.  The  lowest 
frequencies  ( f < . 0 7 )  frequently  show  a  southerly  approach 
direction.  This  may  indicate  the  presence  of  southern  hemisphere 
swell  but  is  also  consistent  with  the  increasing  importance  of 
the  banks  refractive  peak  with  decreasing  frequency.  Note  that 
the  measured  wave  directions  at  the  tower  must  be  modified  by 
refraction  prior  to  comparison  with  SAR  image  collected  in 


deeper  water. 


For  convenience,  Che  values  plotted  in  Figures  8  and  10  for 
6  days  (17,  31  October,  4,  7  Nov,  19,  27  March)  with  Convair  990 
overflights  are  listed  in  Table  3.  'FREQ'  is  the  wave  frequency 
in  hz.  The  column  labeled  "ANG  MAXPWR"  is  the  angle  (relative  to 

270°)  with  the  maximum  energy.  The  angles  listed  are  relative  to 
270°,  i.e.  10°  *  280°,  and  -10°  =  260°.  The  column  labeled  "+- 

deg  to  Halfpwr"  gives  degrees  from  the  peak  angle  to  the  half 
power  points.  The  estimated  spectra  are  usually  quite  symetric. 
The  "FWHM”  is  the  full  width  at  half  maximum  power  in  degrees 
(the  sum  of  the  H —  degrees  column).  The  "power"  listed  is  the 
energy  in  each  frequency  band  and  is  not  related  to  directional 
quantities . 
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Figure  3.  Measured  directional  spectra  at  Torrey  Pines  Beach 

(refracted  into  local  deep  water)  for  10  sample  days 
in  March  1977.  The  data  runs  have  a  cumulative 
length  of  roughly  64  hours  (Pawka,  1983]. 


LOCAL  DEEP  WATER 
MISSION  DAY  ENTRANCE 


Figure  4.  The  theoretical  deepwater  responses  at  Mission  Bay 

Entrance  to  two  north  swell  events.  The  directional 
distributions  of  the  deep  ocean  events  (not  shown) 
have  widths  of  20°  (full  width  at  half  maximum)  and 
are  smooth  (cos  power)  functions.  The  responses  have 
been  block  averaged  into  5°  increments  (Pawka  and 
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Figure  5 . 


The  theoretical  response  at  Mission  Bay  Entrance 
('model  spectrum')  to  a  deep  ocean  mixed  north  swell 
(310°  T)  and  south  swell  (230°  T)  condition.  The 
model  spectrum  (coastal  response)  has  been  blocked 
averaged  into  5°  increments.  Analysis  is  for  a  wave 
period  of  15.0  seconds.  The  smearing  of  the  model 
spectrum  by  finite  length  linear  arrays  is  indicated 
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Figure  6.  Various  maximum  likelihood  estimates  (MLE,  IMLEP , 
IMLEO)  of  a  true  bimodal  spectrum  [Oltman-Shay  and 
Guza ,  1984]. 
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Figure  7.  Estimated  peak  angles  of  bimodal  test  spectra  with 
varying  angular  separation  of  the  modes.  One  mode 
(Ml)  has  half  the  energy  of  the  other  (M2).  Both 
modes  are  10  degrees  full  width  at  half— maximum  power 
[Oltman-Shay  and  Guzn,  1984]. 
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Height  Spectra  Cor  Four  Overflight  Days  in  Fhase-l. 
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Figure  8c-  Wave  Height  Spectra  for  Eight  Overflight  Days  in  Phase -II, 
( Continued) 
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Figure  10a.  Wave  Direction  (at  Maximum  Power)  and  Angular  Spread 
(at  Half  Maximum  Power)  vs.  Frequency,  phase-1. 
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Figure  10b.  Wave  Direction  (at  Maximum  Power)  and  Angular  Spread 
(at  Half  Maximum  Power  Vn.  Frequency),  in  Phase-II. 
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Figure  10c.  Wave  Direction  (al  Maximum  I’ower)  and  Angular  Spread 
(at  Haif  Maximum  Power  Vn.  Frequency).  (Continued) 


Table  1  .  Environmental  Conditions  During  Flights  in  October-November ,  1984: 
Phase-I.  Noted  Times  are  in  Pacific  Daylight  Time  and  Pacific 


Standard  Time  Before  and 

After  28  October  1984,  Respectively. 

Date 

17  Oct  1984 

31  Oct  1984 

4  Nov  1984 

7  Nov  1984 

Time  Analyzed 

1250  -  1413 

1300  -  1600 

1512  -  1625 

1412  -  1625 

Wind  Speed  (knots) 

15 

6-8 

5-7 

10 

Wind  Direction  (deg) 

-275 

-280 

-300 

-250 

Mean  Current  (cm/sec) 
(5m  below  surface) 

<5 

<5 

<5 

-10 

Significant  Wave 
Height  (cm) 

126 

124 

158 

110 

Air  T  ( °C) 

17.2 

17.5 

18.5 

17.5 

1-18 


Vt 

Table  3a.  17 

October 

1984,  1250 

-1413,  40 

degrees 

of  f  reedom . 

. 

■*: 

'v 

PR  PO  AK'G 

MAX PUR 

+ (-)DEO 
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‘■v. 
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0.  7 OF  H  04 

*“  i 

V. 

0. 190E+00 

B.  00 

14 

1  5 

28.  00 

0.  721-H  04 

v| 

0. 194E+00 

2.  00 

23 

25 

47.  00 

0.  5 1  Eh  04 

0. 198E+00 

1.  00 

33 

39 

71 .  00 

0.  4  5E-I04 

+-v^ .» 

0. 202E+00 

8.  00 

15 

1  5 

29.  00 

0.  4  1  Eh  04 

0.  206E+00 

5.  00 

16 

3  7 

32.  00 

0.  63FH04 

0.  209E4  00 

0.  00 

17 

18 

34.  00 

0.  52F.H04 

0.  213E+00 

1  1 . 00 

14 

1  5 

28.  00 

0  47FH04 

0.  217E400 

4 , 00 

19 

21 

39.  00 

0.  39EH04 

'  u 

0.  221 Eh  00 

6.  00 

23 

24 

46.  00 

0.  54 Eh  04 

■«  . 

0. 225E+00 

10.  00 

16 

16 

31. 00 

0.  47EH04 

..\fl 

0. 229E400 

7.  00 

16 

17 

32.  00 

0.  45F.H  04 

0. 233E+00 

19.  00 

18 

19 

36.  00 

0.  37EH04 

0. 237E+00 

0.  00 

20 

e  c. 

41 . 00 

0.  39EH04 

0. 241E+00 

8.  00 

18 

19 

36.  00 

0.  36F.H04 

0. 245E+00 

19.  00 

23 

85 

47.  00 

0.  30EH04 

0. 249E+00 

7.  00 

22 

23 

44.  00 

0.  38EH04 

>v 

0. 252E+00 

12.  00 

32 

33 

64.  00 

0.  24Eh  04 

0. 256E+00 

3.  00 

22 

23 

44.  00 

0.  26Eh  04 

0. 260E+00 

16.  00 

26 

86 

51.  00 

0.  21  Eh  04 

'V 

0. 264E+00 

1  1.  00 

39 

44 

82.  00 

0.  25E  *  01 

0. 26BE+00 

10.  00 

30 

30 

59.  00 

0.  23E  h  04 

0. 272E+00 

9.  00 

23 

25 

47.  00 

0.  34EH04 

0. 276E+00 

9.  00 

31 

32 

62.  00 

0.  25EH04 

0. 280E+00 

2.  00 

36 

38 

73.  00 

0.  21  Eh  04 

0.  28'IE-*  00 

12.  00 

4  5 

47 

93.  00 

0.  22EH04 

0.  28BE  4 00 

12.  00 

35 

37 

7  1 .  00 

0.  23E n  04 

0. 292E+00 

5.  00 

33 

33 

65.  00 

0.  15EH04 

0  295EH00 

-2.  00 

24 

84 

4  7  .  00 

0.  20F  h  04 

O 

3  OO 

ie 

1-20 

20 

37.  00 

0.  1 8F  h  04 

rj  r 

Table  3b.  31  October  1984,  1300-1600,  88  degrees  of  freedom. 


FREQ  ANG  MAXPWR 

+ (“ ) DEG 

0.  493E-01 

-15.  00 

37 

0. 532E-01 

-18.  00 

19 

0. 571E-01 

“25.  00 

13 

0.  610E-01 

“27.  00 

14 

0.  649E-01 

“22.  00 

19 

0. 688E-01 

“10.  00 

12 

0. 728E-01 

-9.  00 

14 

0. 767E-01 

“10.  00 

22 

0. 8C6E-01 

-5.  00 

25 

0. S45E-01 

“15.  00 

23 

0.  884E-01. 

-3.  00 

21 

0.  923E-01 

7.  00 

16 

0. 962E-01 

10.  00 

10 

0.  100E+00 

12.  00 

10 

0.  104E+00 

13.  00 

8 

0.  108E+00 

14.  00 

8 

0.  112E+00 

13.  00 

10 

0.  116E+00 

1 1. 00 

8 

0.  120E+00 

13.  00 

9 

0.  124E+00 

14.  00 

13 

0.  127E+00 

11.  00 

9 

0.  131E+00 

13.  00 

9 

0.  135E+00 

12.  00 

9 

0.  139E+00 

12.  00 

1  1 

0.  143E+00 

14.  00 

12 

0.  147E+00 

12.  00 

13 

0.  151E+00 

13.  00 

13 

0.  155E+00 

10.  00 

12 

0. 159E+00 

8.  00 

13 

0. 163E+00 

10.  00 

14 

0.  167E+00 

12.  00 

13 

0.  1  70E+00 

9.  00 

19 

0.  174E+00 

15.  00 

19 

0.  178E+00 

6.  00 

17 

0.  182E+00 

9.  00 

15 

0.  186E+00 

12.  00 

17 

0.  190E+00 

6.  00 

17 

0.  194E+00 

8.  00 

23 

0.  198E+00 

4.  00 

26 

0.  202E+00 

5.  00 

17 

0. 206E+00 

-1.  00 

17 

0.  209E+00 

4.  00 

24 

0.  213E+00 

4.  00 

23 

0.  217E+00 

-1.  00 

18 

0. 221E+00 

0.  00 

25 

0.  225E+00 

-5.  00 

21 

0.  229E+00 

-3.  00 

20 

0. 233E+00 

-2.  00 

25 

0.  237E+00 

4.  00 

23 

0. 241E+00 

“2.  00 

30 

0.  245E+00 

“1.  00 

20 

0. 249E+00 

-9.  00 

20 

0. 252E+00 

“10.  00 

22 

0.  256E+00 

“1.  00 

25 

0.  260E+00 

-1.  00 

21 

0.  264E+00 

.  -3.  00 

22 

0.  268E+00 

1. 00 

27 

0.  272E+00 

3.  00 

30 

0.  276E+00 

“1 1.  00 

35 

0.  280E+00 

-7.  00 

27 

0. 284E+00 

-7.  00 

22 

0. 2BBE+00 

-6.  00 

34 

0. 292E+00 

-4.  00 

25 

0. 295E+00 

“10.  00 

28 

0.  299E  +  00 

-3.  00 

28 

3  HALFRWR 

FWHM 

POWER 

39 

75.  00 

0.  16E+03 

19 

37.  00 

0.  41E+03 

14 

26.  00 

0. 16E+04 

1  5 

28.  00 

0.  55E+03 

19 

37.  00 

0.  43E+03 

13 

24.  00 

0.  1  1E+04 

14 

27.  00 

0.  19E+04 

21 

42.  00 

0.  30E+04 

22 

46.  00 

0.  20E+04 

25 

47.  00 

0.  17E+04 

20 

40.  00 

0.  93E+03 

1  7 

32.  00 

0.  1  3E+04 

1  1 

20.  00 

0.  30E+04 

1 1 

20.  00 

0.  48E+04 

9 

16.  00 

0.  68E+04 

9 

16.  00 

O.  10E+05 

12 

»  21.  00 

0.  94E+04 

9 

16.  00 

0..1  5E+05 

1  1 

19.  00 

0.  16E+05 

14 

26.  00 

0.  86E+04 

10 

18.  OO 

0.  12E+05 

10 

18.  00 

0.  96E+04 

1  1 

19.  00 

0.  12E+05 

12 

22.  00 

0.  10E+05 

12 

23.  00 

0.  78E+04 

1  4 

26.  00 

0. 88E+04 

14 

26.  00 

0.  86E+04 

1  4 

25.  00 

0.  83E+04 

13 

25.  00 

0.  65E+04 

1 6 

29.  00 

0.  56E+04 

1  5 

27.  00 

0.  53E+04 

19 

37.  00 

0.  44E+04 

21 

39.  00 

0.  40E+04 

18 

34.  00 

0.  42E+04 

17 

31. 00 

0.  50E+04 

19 

35.  00 

0.  41E+04 

18 

34.  00 

0.  33E+04 

24 

46.  00 

0.  21E+04' 

27 

52.  00 

0.  31E+04 

18 

34.  00 

0.  29E+04 

19 

35.  00 

0.  24E+04 

24 

47.  00 

0.  25E+04 

24 

46.  00 

0.  23E-+04 

20 

37.  00 

0.  26E+04 

26 

50.  00 

0.  14E+04 

23 

43.  00 

0.  16E+04 

21 

40.  00 

0.  14E+04 

26 

50.  00 

D.  15E+04 

24 

46.  00 

0.  13E+04 

30 

59.  00 

0.  13E+04 

21 

40.  00 

0.  13E+04 

21 

40.  00 

0.  1  1E+04 

25 

46.  00 

0.  88E+03 

26 

50.  00 

D.  1  1E+04 

22 

42.  00 

0.  86E+03 

24 

45.  00 

0.  B1E+03 

27 

53.  00 

D.  89E+03 

32 

61. 00 

0.  77E+03 

38 

72.  00 

0.  78E+03 

29 

55.  00 

0.  70E+03 

23 

44.  OC 

0.  65E+03 

34 

67:  OC 

0.  B2E+03 

26 

50.  OC 

0.  77E+03 

29 

56.  OC 

0.  83E  +  03 

30 

57.  00 

0.  52E+03 

Table  3c. 


4  Novemeber  1984,  1512-1625,  40  degrees  of  freedom. 


FROG  AK'G 

MAX PUR 

+  ( - >  DEG 

TO  HALFPWR 

FWHI1 

0. 4930-01 

-9.  00 

12 

13 

24.  00 

0.  5320-01 

-5.  00 

9 

10 

18.  00 

0.  5710-01 

2.  00 

7 

8 

1  4.  00 

0. 6100-01 

3.  00 

8 

10 

17.  00 

0. 6490-01 

-2.  00 

10 

1  1 

20.  00 

0. 60GE-0 1 

-6.  00 

1  1 

12 

22.  00 

0.  7280-01 

0.  00 

10 

12 

21 . 00 

0. 7670-0 1 

0.  00 

12 

13 

24.  00 

0  8060-01 

3.  00 

9 

1  1 

19.  00 

0.  84  50-01 

2.  00 

12 

12 

23.  00 

0. e84E-01 

3.  00 

1  2 

12 

23.  00 

0  923E-01 

1 . 00 

13 

13 

25.  00 

0  962E-0 1 

7.  00 

B 

10 

17.  00 

0.  1000  +  00 

4.  00 

10 

1  1 

20.  00 

0. 1040+00 

B.  00 

8 

10 

17.  00 

0.  106E+00 

6.  00 

9 

10 

18.  00 

0.  1  12E+00 

13.  00 

8 

10 

17.  00 

0  116E+00 

7.  00 

10 

10 

19.  00 

0.  120E  +  00 

7.  00 

9 

1  1 

19.  00 

0  124E+00 

7.  00 

1  1 

12 

22.  00 

0  1270+00 

1  1  .  00 

12 

1  2 

23.  00 

0. 1310+00 

8.  00 

10 

1  1 

20.  00 

0.  1350  +  00 

6.  00 

1  6 

16 

31 . 00 

0. 139E+00 

3.  00 

14 

1  4 

27.  00 

0. 143E+00 

13  00 

1  5 

1  5 

29.  00 

0.  147E+00 

10.  00 

13 

14 

26.  00 

0.  151E+00 

7.  00 

21 

22 

42.  00 

0  155E+00 

- 1 . 00 

12 

14 

25  00 

0. 1590+00 

10  00 

1  1 

12 

22.  00 

0. 1630+00 

4.  00 

1  5 

16 

30.  00 

0.  1670+00 

4.  00 

1  5 

17 

31 . 00 

0.  1700+00 

1 . 00 

14 

14 

27.  00 

0.  1740  +  00 

9.  00 

21 

23 

43.  00 

0. 1780+00 

0.  00 

17 

ie 

34.  00 

0  1G2E+00 

6.  00 

20 

21 

40.  00 

0.  1860*00 

-2.  00 

17 

18 

34.  00 

0.  1900+00 

6.  00 

16 

1  7 

32.  00 

0.  1940  +  00 

0.  00 

16 

1  7 

32.  00 

0.  1980  +  00 

-1.  00 

16 

1  7 

32.  00 

0. 2020+00 

1  1. 00 

21 

21 

4  1 . 00 

0. 206E+00 

1.  00 

23 

23 

45.  00 

0  2090*00 

8  00 

33 

34 

66.  00 

0  2130*00 

0.  00 

21 

O  ^ 
c.  e. 

42.  00 

0.  21 70*  00 

0  00 

1  7 

ie 

34.  00 

0.  2210  +  00 

6.  00 

1  5 

1  7 

31 . 00 

0  2250*00 

-4.  00 

17 

19 

35.  00 

0  2290+00 

2.  00 

1  7 

18 

34.  00 

0.  2330  +  00 

-2.  00 

1  6 

17 

32.  00 

0.  237E+00 

6.  00 

21 

23 

43.  00 

0  241E+00 

-2.  00 

21 

23 

43.  00 

0  245E+00 

2.  00 

18 

20 

37.  00 

0  249E+00 

5.  00 

20 

21 

40.  OC 

O'.  252E+00 

0.  256E+00 

-1. 00 

16 

17 

32.  00 
****** 

0.  260E+00 

0.  00 

17 

19 

35.  00 

0.  264E+00 

3.  00 

21 

21 

4  1 . 00 

0.  268E+00 

0.  272E+00 

0.  276E  +  00 

0.  280E  +  00 

0.  204E+00 

0. 20GE+00 

0. 292E+00 

4.  00 

1  7 

19 

35.  00 
****** 
****** 
****** 
****** 

****** 

****** 

0. 2950+00 

0.  2990  +  00 

1 . 00 

1  3 

14 

26.  00 

*«*««« 

POWER 
0.  420*04 
0.  1  5E  *  Op 
O  54F.+  05 
O.  DIE*  00 
0.  230*05 
0.  190*05 
O.  210*05 
O.  1  4E  •  05 
O  110*05 
O.  120*05 
O.  000*04 
O.  060*04 
O.  73E  +  04 
O.  75E  +  04 
0.  1  4E*05 
0.  1  10+05 
0.  1  1E  +  05 
0.  1  10*05 
0.  1  10*05 
0.  360*04 
O.  600*04 
0  420*04 
O.  500*04 
O.  34  0*04 
0.  200*04 
0.  300*04 
0  270+04 
0.  340  +  04 
0.  370*04 
0.  290  +  04 
0.  1  90  +  04 
O.  290*04 
0.  210*04 
0.  100*04 
0.  130  +  04 
0.  130  +  04 
0.  170  +  04 
0.  190*04 
0.  14  0*04 
0.  140*04 
0.  990+03 
O.  120  +  04 
0.  1  10*04 
O.  120+04 
0.  1  1  0*  04 
0.  130+04 
0.  130  +  04 
0.  130+04 
0.  100*04 
0.  130*04 
0.  930  +  03 
0.  620*03 
0.  140+04 
0.  760*03 
0.  13E+04 
0.  120*04 
0.  190  +  04 
0.  140*04 
0.  1  70*04 
0.  1  70*  04 
0.  200+04 
0.  210+04 
0.  260*04 
0.  200*04 
O.  300 +  04 


Table  3d.  7  November  1984,  1412-1625,  40  degrees  of  freedom. 


FREQ  ANG 

HAXFTs’R 

0. 493E-01 

-18.  00 

0.  532E-0 1 

-11.  00 

0.  571E-01 

-13.  00 

0. 61CE-01 

-15.  00 

0. 649E-01 

-20.  00 

0. 688E-01 

-8.  00 

0. 72BE-01 

1 .  00 

0. 767E-01 

5.  00 

0. 806E-01 

8.  00 

0. B45E-01 

6.  00 

0.  8 8 4 E - 0 1 

6.  00 

0. 923E-01 

7.  00 

0. 962E-01 

5.  00 

0.  lOC’E+OO 

10.  00 

0. 104E+00 

4.  00 

o.  iceE-too 

7.  00 

0.  112E-+00 

4.  00 

0. 1 16E+00 

9.  00 

0.  120E  +  00 

4.  00 

0. 124E+00 

7.  00 

0.  1 2  7E-<  00 

4.  00 

0.  1  3 1 E  <00 

8.  CO 

0.  1  35E  *  00 

8.  00 

0  !39E^0 

19.  00 

0.  143EHOO 

16.  00 

0.  1  4  7E-*  00 

13.  00 

0.  151E400 

22.  00 

0.  155E^0 

20.  00 

0.  1  59E-+  00 

20.  00 

0.  163L  +  00 

26.  00 

0.  1 67E  <  00 

22.  00 

0.  1 70E '  00 

26  00 

0.  1  74E  <  00 

18.  00 

0.  1  76E-<  00 

21. 00 

0.  182E+00 

23.  00 

0.  186E^0 

24.  00 

0.  1 90E-<  00 

15.  00 

0.  1 94E  <  00 

13.  00 

0.  1 9BE<-00 

19.  00 

0.  202E-<  00 

11  .  00 

0.  206EH0O 

1  1 . 00 

0.  2C9E-<  00 

5.  00 

0.  21  3F  <  00 

2.  00 

0. 217E  <00 

18.  00 

0.  22 1  E<  00 

0.  00 

0.  225E-*  00 

9.  00 

0.  229E-*  00 

12.  00 

0. 233E+00 

14.  00 

0.  237E  +  00 

4.  00 

0. 241E+00 

8.  00 

0.  245E-<00 

2.  00 

0  249E+00 

6.  00 

0.  252E-<  00 

2.  00 

0.  256E  +  00 

13.  00 

0.  260E  +  00 

7.  00 

0.  264 E  <  00 

10.  00 

0.  268E+00 

6.  00 

0.  272E-<  00 

1. 00 

0.  276E->  00 

-6.  00 

0.  280E-<00 

-2.  00 

0.  284E+00 

3.  00 

0.  28SE-<  00 

-2.  00 

0.  292E  <00 

-14.  00 

0.  295E-<  00 

-1 1 . 00 

0  295E-<00 

-1 1 . 00 

+  (  -  )  DEG  10  MALFPIJR 


17 

16 

22 

'"I 

c  c. 

13 

13 

10 

1  1 

15 

19 

18 

1  7 

10 

1  2 

9 

12 

10 

1  1 

1  1 

13 

9 

1  0 

10 

1  0 

12 

1  4 

9 

1  1 

12 

13 

18 

25 

14 

1  6 

13 

1  4 

14 

1  4 

15 

1  7 

17 

on 
c  c. 

16 

38 

16 

1  7 

16 

3  9 

12 

14 

15 

1  7 

1  1 

1  2 

12 

1  4 

1  1 

12 

12 

12 

1  1 

13 

13 

1  4 

18 

"n 

e_  r. 

24 

26 

19 

23 

19 

21 

23 

26 

24 

8  / 

25 

29 

34 

4  4 

34 

-10 

30 

30 

35 

36 

4  2 

6  7 

43 

39 

39 

42 

36 

40 

25 

25 

22 

rn 

C-  c. 

2? 

31 

30 

31 

29 

31 

27 

2  7 

40 

42 

32 

33 

43 

4  5 

30 

31 

49 

56 

41 

45 

33 

33 

57 

54 

4  5 

47 

55 

65 

59 

61 

56 

5  7 

FWIIM  POWER 

32.  00  0.  2  1  E-»  03 
'13.  00  O.  3  2FHOJ 

25.  00  ■  0.  51E  +  03 

20.  00  0.  29E+04 

33.  00  •  0.  12E+05 

34.  00  0.  1  3E-»05 

21.  00  0.  3  6EH05 

20.  00  0.  1  05 

20.  00  0.  1  4E  +  05 
23.00  0.  86E-»04 
ie.  00  0.  9  4  E  -<  0  4 
19.00  0.  46E  *04 

25.  00  0.  32E-*  04 
19.  00  0.  49E-*  04 

24.  00  0.  23E+04 

42.  00  O.  22E-<  04 
29.  00  0.  21E-*  04 

26.  00  0  1  8E-*  04 

27.  00  O.  21E  <04 
31 . 00  O.  16E-<04 
3B.  00  O.  1  2Eh  04 

33.  00  0.  1  IF. -i  04 
32.  00  0.  1  4E  <  04 

34.  00  0.  1  2E  <  04 

25.  00  0.  1  8E-<  04 
31 . 00  0.  10En>4 

22.  00  0.  24E^04 

25.  00  0  26E-*  04 

22.  00  O.  34EM>4 

23.  00  O.  34E-*  04 
23.  00  0.  4 7 F. -*  0 4 

26.  00  0.  27E -<  04 
39.  00  0.  32E -<04 

49.  00  0.  ieE  +  04 
41.  00  0.  33E-<  04 
39.  00,0.  37E  +  04 

48.  00  0.  261' -<  04 

50.  00.0.  30 E  <  04 
53  00  0.  21E  <04 
77.  00.0.  19E-I04 
73.  00  0.  20F.-<04 
59.  00  0.  22E  <  04 
70.  00  0.  25F.  <  04 

108.  00  O  1  41  i  04 
81. 00,0.  13EU>4 
80.  00|0.  1  2E-<  04 
75.  00,0.  1  5E-<  04 

49.  00  0.  1  3E  <  04 

43.  00,0.  13E-<04 

59.  00,0.  1  6E-<  04 

60.  00(0.  98E-*  03 

59.  00  O.  1  4E-*  0  < 
53.  00,0.  1  3E  ‘ '. 

81 . 00  0.  JOE-1  04 

64.  00  0.  3  2E-*  04 
87.  00|0.  95E-*  03 

60.  00,0.  1  lE-<  04 
104.  00,0.  93E-*  03 

85.  00  0.  68E-*  03 

65.  00;o.  1  0E-*  04 
110.  00  0.  89 E-*  03 

91. 00.0.  6 4 F -t  0 3 
1  19.  00,0.  961. -*03 
1  19.  00  0  7  IF-*  03 
1  12.  00  O.  6PE-»03 
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EREQ  ANG 

riAxrwn 

+ ( — ) DEG 

It)  IIALPPWR 

EMI  HI 

0. 973E-01 

-6.  00 

7 

7 

17. 

00 

0. 532E-01 

-19.  00 

ID 

1  7 

31. 

00 

0. 571E-01 

-18.  GO 

12 

13 

24. 

00 

0. 610E-01 

-14.  00 

1  1 

1  1 

21. 

00 

0.  697E-01 

-17.  00 

13 

14 

26. 

00 

0.  688E-0 1 

-12.  00 

8 

e 

15. 

00 

0. 728E-01 

-4.  GO 

ID 

19 

2B. 

00 

0.  767E-01 

0.  00 

B 

10 

17. 

00 

0. 806E-01 

1. 00 

7 

7 

1  5. 

00 

0.  89  5E-01 

2.  00 

7 

7 

1  5. 

00 

0.  BB9E-01 

2.  00 

B 

B 

1  D. 

00 

0.  723E-01 

3.  00 

1  1 

13 

23. 

00 

0.  762E-01 

3.  00 

? 

10 

18. 

00 

0.  100E+00 

4.  00 

10 

10 

17. 

00 

0.  104E+00 

4.  00 

7 

1  1 

17. 

00 

0. 1C8E+00 

5.  00 

1  1 

13 

23. 

00 

0.  1  12E+00 

5.  00 

r 

7 

1  1 

17. 

00 

0.  1 1  6E+00 

5.  00 

10 

1  1 

20. 

00 

0.  120E+00 

5.  00 

1  0 

1  2 

21  . 

00 

0.  124E+00 

B.  00 

1 1 

13 

23. 

00 

0.  127E+00 

7.  00 

1 1 

13 

23. 

00 

0.  131E+00 

1  1.  00 

7 

10 

18. 

00 

0.  135E+00 

7.  00 

8 

7 

16. 

00 

0.  137E+00 

7.  00 

10 

1  1 

20. 

00 

0.  193E+00 

12.  00 

12 

13 

24. 

00 

0.  147E+00 

7.  00 

12 

13 

24. 

00 

0.  151E+00 

7.  00 

10 

13 

22. 

00 

0.  155E+00 

15.  00 

3  3 

15 

27. 

00 

0.  157E+00 

4.  00 

1  1 

12 

22. 

00 

0.  163E+00 

7.  00 

18 

18 

35. 

00 

0.  1  67E+00 

7.  00 

16 

10 

33. 

00 

0.  1  70E+00 

6.  00 

17 

18 

34. 

00 

0.  174E+00 

6.  00 

19 

16 

27. 

00 

0.  178E+00 

6.  00 

13 

1  4 

26. 

00 

0.  182E+00 

3.  00 

17 

18 

34. 

00 

0.  1  BfeE+OO 

11.00 

16 

1  7 

32. 

00 

0.  170E+00 

1.  00 

16 

16 

31  . 

00 

0.  1 74E^  00 

4.  00 

20 

22 

41. 

00 

0.  178E+00 

3.  00 

17 

O 

C-  *=. 

40. 

00 

0.  202E+00 

5.  00 

22 

24 

45. 

00 

0. 206E+00 

6.  00 

16 

17 

34. 

00 

0.  207E+00 

7.  00 

17 

17 

35. 

00 

0.  213E+00 

B.  00 

17 

20 

38. 

00 

0.  217E+00 

2.  00 

25 

e-  » 

53. 

00 

0.  221E+00 

3.  00 

21 

e-c. 

42. 

00 

0.  225E+  00 

-13.  00 

32 

35 

66. 

00 

0. 227E+00 

6.  00 

31 

36 

66. 

00 

0. 233E+00 

1  1. 00 

17 

21 

37. 

00 

0. 237E+00 

4.  CO 

31 

33 

63. 

00 

0.  241E+00 

4.  00 

22 

24 

45 

00 

0. 245E+00 

- 1 . 00 

25 

26 

50. 

00 

0.  247E-*  00 

2.  00 

20 

22 

41 . 

00 

0.  252E  +  00 

2.  00 

32 

32 

63. 

00 

0.  256E+00 

4.  00 

26 

27 

52. 

00 

0.  260E+00 

5.  00 

27 

27 

53. 

00 

0.  264E-*00 

0.  00 

30 

31 

60. 

00 

0.  26BE+00 

B.  00 

27 

30 

58. 

00 

0.  272E+00 

-3.  00 

3D 

37 

71 . 

00 

0. 276E400 

-4.  00 

33 

34 

66. 

00 

0.  2B0E-*  00 

1. 00 

30 

31 

60. 

00 

0.  2B9E-*00 

4.  00 

31 

33 

63 

00 

0.  2BBE  +  00 

6.  00 

39 

39 

67 

00 

0.  272E*  00 

B.  00 

27 

31 

57. 

00 

0.  275E-»00 

3.  00 

37 

37 

75. 

00 

O.  277E*  00 

-5.  OO 

31 

32 

62. 

00 

power 
o.  1  OF. 1  09 

o.  2?e •«  oj 

0.  7  1 E-*  03 
0.  62E+03 

o.  ue+04 
0.  67E^  09 
0.  84E  *04 
0.  16E+05 
0.  15E-*05 
0.  13E.-*05 
0  12E^05 
0.  62E  *  09 
0.  7 1 E  *  04 
0.  B2E+04 
0.  7  1  E  ■*  0  4 
0.  D?E»04 
0.  0?E  *  09 
0.  7  IE  *09 
0.  ODE  '  09 
0  78E  *  09 
0.  8  7E-*  09 
0.  1  4E-*  05 
0  1PE'05 
0  1  6E  *  05 
0  1  1E^  05 
0  7 7E -*  09 
O  1 OE  *  05 
0  7  5E *  09 
0.  76E-*  09 
0  5  IE* 09 
0  6 9 E •'  09 
0  40E*04 
O.  4  6E*  09 
0  45E^  09 
0.  'JOE  *09 
0  9  3E-*  09 
0.  4  7  E  ■*  0  4 
0.  27E  *  09 
0.  9  4E-*  09 
0.  35E^04 
0.  32EH  09 
0.  20E  ♦  09 
O.  23E-*  04 
0.  31E-*04 
0.  1BE  +  09 
0.  17E+04 
0. 23E+04 
0.  24E+04 
0.  19E  +  09 
0  17E+04 
0.  14E+04 
0.  15E+09 
0.  1  0E-*  09 
0.  1  5E^  09 
0-  12E+04 
0.  1  1E*04 
0.  B7E  +  03 
0.  12E+04 
0.  BBE+03 
0.  10E+09 
0.  88E+03 
0.  80E+03 
V  1  1E  +  04 
0  73E-*03 
O.  50E+O3 
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FREQ  ANG 

NAXPWR 

+1-1UEG 

0. 493E-01 

-7.  00 

27 

0.  532E-01 

-11.  00 

18 

0. 571E-01 

-IB.  00 

1  5 

0.  610E--01 

-8.  00 

1 6 

0. 649E-01 

-8.  00 

12 

0. 68BE-01 

1.  00 

14 

0. 72EE-01 

-2.  00 

12 

0. 767E-01 

-3.  00 

12 

0. B06E-01 

-4.  00 

12 

0.  S45E-01 

1.  00 

12 

0.  8S4E-01 

3.  00 

8 

0. 923E-01 

3.  00 

10 

0. 962E-01 

4.  00 

10 

0.  100E+00 

3.  00 

9 

0.  104E+00 

4.  00 

10 

0. 108E+00 

9.  00 

9 

0.  112E+00 

6.  00 

9 

0. 116E+00 

8.  00 

10 

0. 120E+00 

8.  00 

1  1 

0.  124E+00 

8.  00 

13 

0.  127E+00 

11. 00 

9 

0. 131E+00 

6.  00 

10 

0.  135E+00 

8:  00 

12 

0. 139E+00 

13.  00 

1  4 

0.  143E+00 

11.  00 

14 

0.  147E+00 

11.  00 

1 6 

0. 151E+00 

8.  00 

1  6 

0.  155E+00 

12.  00 

14 

0.  159E+00 

4.  00 

17 

0.  163E+00 

7.  00 

23 

0.  167E+00 

6.  00 

1 6 

0.  170E+00 

13.  00 

17 

0.  174E+00 

13.  00 

17 

0. 17BE+00 

B.  00 

24 

0. 1B2E+00 

7.  00 

1  6 

0. 186E+00 

15.  00 

17 

0.  190E+00 

3.  00 

17 

0.  194E+00 

5.  00 

1  5 

0.  19GE+00 

8.  00 

18 

0.  202E+00 

6.  00 

1  7 

0.  206E+00 

5.  00 

19 

0.  2C9E+00 

3.  00 

20 

0.  213E+00 

1.  00 

19 

0.  217E+00 

3.  00 

22 

0.  221E+00 

8.  00 

25 

0.  225E+00 

-1. 00 

18 

0.  229E+00 

-6.  00 

25 

0. 233E+00 

-3.  00 

25 

0.  237E+00 

0.  00 

23 

0. 24 1E+00 

-1.  00 

22 

0.  245E+00 

-16.  00 

26 

0. 249E+00 

2.  00 

30 

0.  252E+00 

-1.  00 

24 

0.  256E+00 

-7.  00 

30 

0. 260E+00 

-11. 00 

45 

0. 264E+00 

- 1 1 .  00 

29 

0.  268E+00 

-5.  00 

30 

0. 272E+00 

-15.  00 

44 

0. 276E+00 

-16.  00 

46 

0.  280E+00 

-6.  00 

37 

0. 284E+00 

1. 00 

39 

0.  288E+00 

-1.  00 

4  1 

0.  292E+00 

-8.  00 

39 

0. 295E+00 

-12. 00 

4  1 

0. 299E+00 

-9.  00 

35 

MAl.FPWR 

rwnu 

PUt  HR 

27 

53.  00 

0.  4  31  '02 

18 

35.  00 

0.  6't  '02 

15 

29.  00 

0.  1  it:«03 

17 

32.  00 

0.  16E4  03 

13 

24.  00 

0.  51E403 

15 

28.  00 

0.  24F.+  04 

14 

25.  00 

0.  41E+04 

12 

23.  00 

0.  65E+04 

13 

24.  00 

0.  62E '  04 

13 

24.  00 

0.  35E  *  04 

9 

16.  00 

0.  55E+04 

10 

19.  00 

0.  29E  +  04 

12 

21. 00 

0.  31 E+04 

10 

18.  00 

0.  29F.4  04 

10 

19.  00 

0.  30E4  04 

1  1 

19.  00 

0.  30E  +  04 

10 

18.  00 

0.  30F.-'  04 

10 

19.  00 

0.  32E+04 

1  1 

21.  00 

0.  25E+04 

1  5 

27.  00 

0.  19E+01 

10 

18.  00 

0.  23E +04 

1  1 

20.  00 

0.  20E+04 

13 

24.  00 

0.  18E+04 

1  5 

2B.  00 

0.  15E+04 

16 

29.  00 

0.  ieE+04 

17 

32.  00 

0.  16E+04 

18 

33.  00 

0.  14E+04 

17 

30.  00 

0.  22E+04 

ie 

34.  00 

0.  16E  +  04 

28 

50.  00 

0.  16E+04 

17 

32.  00 

0.  2BE+04 

19 

35.  00 

0.  27E+04 

20 

36.  00 

0.  36E+04 

32 

55.  00 

0.  31E+04 

16 

31. 00 

0.  40E+04 

18 

34.  00 

0-  51E+04 

19 

35.  00 

0.  44E+04 

16 

30.  00 

0.  60E+04 

21 

38.  00 

0.  50E+04 

18 

34.  OO 

0-  57E.4  04 

21 

39.  00 

0.  61E404 

22 

41.  00 

0.  40E+04 

19 

37.  00 

0.  50E4  04 

24 

45.  00 

0.  40E4  04 

30 

54.  00 

0.  36F.4  04 

20 

37.  00 

0.  31F.4  04 

26 

50.  00 

0.  38E+04 

25 

49.  00 

0.  35E+04 

24 

46.  00 

0.  21E+04 

23 

44.  00 

0.  23E+04 

27 

52.  00 

0.  24E+04 

32 

61.  00 

0.  23E+04 

26 

49.  00 

0.  16E+04 

31 

60.  00 

0.  12E+04 

43 

87.  00 

0.  19E+04 

31 

59.  00 

0.  14E+04 

32 

61. 00 

0.  17E404 

48 

91. 00 

0.  12E+04 

4? 

94.  00 

0.  16E+04 

38 

74.  00 

0.  13E+04 

41 

79.  00 

0.  15E4  04 

40 

BO.  00 

0.  16E+04 

41 

79.  OO 

0.  13C+04 

46 

86.  00 

0.  15E+04 

36 

70.  OO 

0.  1  1F.  +  04 
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Development  and  Utilization  of  a 
Surface  Energy  Measurement  System 
in  TOWARD 

William  D.  Garrett 
Naval  Research  Laboratory 
Washington,  D.  C.  20375 

BACKGROUND: 

Polar  organic  molecules  absorb  at  the  air-sea  interface 
where  they  may  form  continuous  films  which  modify  numerous  small- 
scale  interfacial  processes.  When  the  surface  concentration  of  the 
surface-active  material  attains  a  level  where  the  molecules  of  the 
surface  film  come  into  contact  (about  1  mg  per  sq.  meter),  the 
surface  becomes  relatively  incompressible,  a  thin  layer  of  water  at 
the  air-water  interface  is  immobilized,  and  capillary  and  short 
gravity  waves  are  attenuated.  Consequently,  backscattered  micro- 
wave  radar  signals  are  affected  by  these  film-induced  modifications 
of  the  sea  surface,  a  condition  which  occurs  when  the  surface  tension 
has  been  reduced  by  as  little  as  1-2  mN  per  meter. 

The  effects  of  organic  films  on  small-scale  dynamics  of  the 
air-sea  interface  are  influenced  by  the  film  compressibility  and 
film  flow  across  a  surface  tension  gradient  (Marangoni  effect). 

Surface  tension  is  inherent  in  both  of  these  parameters,  and  is  itself 
an  indicator  of  the  existence  of  a  coherent  organic  film  capable  of 
modifying  interfacial  processes.  For  these  reasons,  it  was  considered 
important  to  include  the  measurement  and  monitoring  of  surface  tension 
in  the  Tower  Ocean  Wave  and  Radar  Dependence  (TOWARD)  Experiment. 

SURFACE  ENERGY  MEASUREMENT  SYSTEM 

A  computer-operated  system  for  the  semi-continuous  measure¬ 
ment  of  the  surface  tension  of  seawater  was  developed,  constructed, 
mounted  on  the  NOSC  Tower,  and  operated  during  the  two  phases  of  the 
TOWARD  Experiment.  The  instrumentation  included  (1)  a  cetrifugal 
pump  mounted  aboard  a  tethered  raft  (fig.  1),  (2)  a  tower-mounted 

surface  tension  (fig.  2),  and  (3)  a  microcomputer.  The  sensor  as 
well  as  the  sampling  pump  were  interfaced  to  the  microcomputer.  At 


predetermined  intervals  seawater  was  pumped  from  the  upper  2  cm  of 
the  ocean  to  the  measurement  chamber.  Surface  tension  was  measured 
automatically  by  a  compu t e r- con t r o 1 1 e d  strain-gauge  sensor  attached 
to  a  hydrophilic  plate  in  contact  with  the  water  in  the  measurement 
chamber.  A  graph  of  the  surface  tension  vs.  time  was  displayed  on 
the  computer's  video  screen  and  updated  every  three  minutes  to 
present  a  real-time  record  of  the  relative  concentration  of 
su r f ac e - a c t i v e  organic  material  in  the  surface  water  passing  the 
tower.  At  the  end  of  each  hour,  a  printer  produced  a  hard  copy 
of  the  graph  of  surface  tension  vs.  time  during  the  previous 
one-hour  period. 

Microcomputer  and  Interfaces:  An  Apple  II-E  microcomputer 
with  64  k  of  RAM  and  one  disk  drive  was  used  to  control  the 
sampling  pump  and  to  take  data  from  the  strain  gauge.  A  computer- 
switchable  5-volt  signal  was  obtained  from  the  game  controller 
IC  socket  of  the  computer  to  operate  an  optically  isolated  relay 
which  controlled  the  AC  power  to  the  pump.  Two  analog-to-digital 
converter  cards  were  used  in  the  microcomputer,  one  to  read  the 
strain  gauge  signals,  and  the  other  contained  a  clock  to  determine 
the  actual  time  of  sampling.  A  10.2-volt  signal  was  taken  from 
the  A/D  strain  gauge  card  to  energize  the  Wheatstone  bridge 
circuit  of  the  strain  gauge. 

Measurement  Strategy:  Figure  3  illustrates  the  typical 
sequence  of  events  during  a  3-minute  measurement  cycle.  The  upper 
part  of  the  figure  depicts  the  voltage  reading  on  the  sensor;  the 
lower  portion  shows  the  corresponding  water  level  in  the  measure¬ 
ment  chamber.  Beginning  with  a  nearly  empty  sample  chamber  to 
which  the  pump  has  just  begun  to  send  water,  a  constant  voltage 
reading  is  observed  due  to  the  plate  hanging  in  air.  The  water 
rises  until  it  contacts  the  plate,  at  which  time  the  plate  is 
pulled  downward  by  the  surface  tension  of  the  seawater  causing  a 
rapid  increase  in  voltage  output.  The  water  continues  to  rise  in 
the  chamber  to  the  level  of  the  overflow  port.  This  results  in  a 
reduction  of  the  voltage  output  of  the  sensor  because  of  buoyancy 
of  the  plate  in  the  seawater.  As  the  water  flushes  through  the 
chamber  and  overflows,  a  noisy,  constant  voltage  signal  results. 
When  the  pump  is  switched  off  by  the  computer,  the  water  drains 
slowly  from  the  chamber,  and  the  buoyancy  effect  decreases  as  less 
of  the  plate  is  submerged.  Once  the  bottom  of  the  sensor  plate  is 


predetermined  intervals  seawater  was  pumped  from  the  upper  2  cm  of 
the  ocean  to  the  measurement  chamber.  Surface  tension  was  measured 
automatically  by  a  computer-controlled  strain-gauge  sensor  attached 
to  a  hydrophilic  plate  in  contact  with  the  water  in  the  measurement 
chamber.  A  graph  of  the  surface  tension  vs.  time  was  displayed  on 
the  computer's  video  screen  and  updated  every  three  minutes  to 
present  a  real-time  record  of  the  relative  concentration  of 
surface-active  organic  material  in  the  surface  water  passing  the 
tower.  At  the  end  of  each  hour,  a  printer  produced  a  hard  copy 
of  the  graph  of  surface  tension  vs.  time  during  the  previous 
one-hour  period. 

Microcomputer  and  Interfaces:  An  Apple  II-E  microcomputer 
with  64  k  of  RAM  and  one  disk  drive  was  used  to  control  the 
sampling  pump  and  to  take  data  from  the  strain  gauge.  A  computer- 
switchable  5-volt  signal  was  obtained  from  the  game  controller 
IC  socket  of  the  computer  to  operate  an  optically  isolated  relay 
which  controlled  the  AC  power  to  the  pump.  Two  analog-to-digital 
converter  cards  were  used  in  the  microcomputer,  one  to  read  the 
strain  gauge  signals,  and  the  other  contained  a  clock  to  determine 
the  actual  time  of  sampling.  A  10. 2-volt  signal  was  taken  from 
the  A/D  strain  gauge  card  to  energize  the  Wheatstone  bridge 
circuit  of  the  strain  gauge. 

Measurement  Strategy:  Figure  3  illustrates  the  typical 

sequence  of  events  during  a  3-minute  measurement  cycle.  The  upper 
part  of  the  figure  depicts  the  voltage  reading  on  the  sensor;  the 
lower  portion  shows  the  corresponding  water  level  in  the  measure¬ 
ment  chamber.  Beginning  with  a  nearly  empty  sample  chamber  to 
which  the  pump  has  just  begun  to  send  water,  a  constant  voltage 
reading  is  observed  due  to  the  plate  hanging  in  air.  The  water 
rises  until  it  contacts  the  plate,  at  which  time  the  plate  is 
pulled  downward  by  the  surface  tension  of  the  seawater  causing  a 
rapid  increase  in  voltage  output.  The  water  continues  to  rise  in 
the  chamber  to  the  level  of  the  overflow  port.  This  results  in  a 
reduction  of  the  voltage  output  of  the  sensor  because  of  buoyancy 
of  the  plate  in  the  seawater.  As  the  water  flushes  through  the 
chamber  and  overflows,  a  noisy,  constant  voltage  signal  results. 
When  the  pump  is  switched  off  by  the  computer,  the  water  drains 
slowly  from  the  chamber,  and  the  buoyancy  effect  decreases  as  less 
of  the  plate  is  submerged.  Once  the  bottom  of  the  sensor  plate  is 
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level  with  the  water  surface,  only  the  surface  tension  force  is 
operative.  As  the  water  level  continues  to  fall,  the  plate  is 
detached  from  the  water  surface,  and  the  maximum  force  is  recorded. 
Subsequently,  a  steauy  voltage  due  to  the  weight  of  the  plate  in 
air  results  and  is  recorded.  The  difference  in  voltage  between 
that  due  to  the  maximum  pull  on  the  plate  and  that  due  to  the 
weight  of  the  plate  in  air  is  directly  proportional  to  the  surface 
tension  of  the  water  in  the  chamber.  Appendix  I  is  a  flow  chart 
of  the  computer  program  which  controls  the  Surface  Energy 
Measurement  System. 

IN  SITU  SURFACE  TENSION  MEASUREMENT  AND  INSTRUMENT  CALIBRATION 

A  series  of  laboratory  calibrated  spreading  oils  were  used 
periodically  to  determine  the  seawater  surface  tension  in  situ  for 
comparison  with  the  Surface  Energy  Measurement  System  (SEMS). 
Special  events  such  as  natural  slicks  associated  with  kelp  or 
internal  waves  were  given  particular  attention.  In  addition,  a 
du  Nouy,  ring-type,  interfacial  tensiometer  was  employed  to 
measure  the  surface  tension  of  collected,  film-free  surface 
seawater  at  constant  temperature  (T)  and  known  salinity  (S),  as 
a  test  of  the  empirical  expression. 

Surface  Tension  =  75.64  -  0.144T  +  0.02215S, 
where  T  is  in  degrees  C  and  S  is  in  ppt.  This  expression  was 
found  to  be  accurate  to  within  +/-  0.1  mN  per  m,  the  accuracy  of 
the  tensiometer.  Thus,  data  on  sea  surface  T  and  S  can  be  used 
to  fine  tune  the  surface  tension  data  determined  by  the  other 
methodologies  employed  in  this  experiment. 

The  SEMS  was  calibrated  and  its  efficacy  verified  by  spread¬ 
ing  know  monomo 1 e c u la r  films  upwind  of  the  floating  sea  surface 
sampling  pump,  so  that  they  would  drift  across  the  unit  and  be 
collected  and  sent  to  the  measurement  chamber.  Figures  4  and  5 
are  computer-generated  plots  of  surface  tension  vs.  time  for 
periods  during  which  the  calibration  slicks  were  deployed.  Surface 
tension  reductions  recorded  agreed  well  with  the  known  values  for 
the  monolayers  applied  to  the  sea  surface.  In  the  case  of  the 
strongly  surface-active  material,  isostearyl  alcohol  (20 E),  the 
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surface  film  was  just  large  enough  to  pass  the  collection  unit 
for  one  pumping  cycle.  There  was  a  residual  effect  of  the  collected 
surface-active  material  for  only  one  subsequent  measurement  cycle, 
after  which  the  system  appeared  to  have  flushed  completely,  and 
the  surface  tension  returned  to  a  value  near  that  of  clean  seawater. 

In  the  case  of  the  glycerol  trileate,  a  very  large  slick  was 
produced,  and  the  effect  of  the  monolayer  was  measured  over  a 
number  of  measurement  cycles.  Two  examples  of  documented  natural 
surface  events  measured  by  SEMS  are  depicted  in  Figures  6  and 
7,  representing  natural  slicks  from  kelp  and  possibly  from  an 
internal  wave  convergence  zone,  respectively.  In  the  case  of  the 
nonkelp  slick  (figure  7),  calibrated  spreading  oil  measurements 
yielded  a  surface  tension  reduction  of  12mN  per  m  at  10-08-1435, 
a  value  close  to  that  measured  by  the  SEMS.  These  results  indicate 
that  the  SEMS  accurately  measures  sea  surface  film  tension  if  the 
sur f ac e-ac t ive  layers  are  sufficiently  large  to  cover  the  sampling 
unit  during  the  pumping  cycle,  even  though  the  system  does  not 
provide  a  direct  in-situ  measurement.  Futhermore,  it  has  been 
demonstrated  that  a  one-molecule-thick  film  can  be  sample  from  the 
sea  surface  and  its  surface  tension  value  accurately  determined  by 
the  SEMS. 

DATA  QUALITY  AND  PLANNED  INTERCOMPARISONS  WITH  OTHER  TOWARD  PRODUCTS 

A  catalog  of  surface  tension  data  taken  with  the  SEMS  for  both 
TOWARD  exper iemental  periods,  20  September  to  3  November  1984,  and 
8-26  April  1985,  is  attached  as  Appendix  II.  These  data  have  been 
normalized  on  the  basis  of  corresponding  in-situ  surface  tension 
measurements  made  with  calibrated  spreading  oils.  The  SEMS  data  will 
be  additionally  fine  tuned  to  account  for  seawater  surface  temperature. 
Although  we  have  temperature  data  for  many  days  of  operation,  it  will 
be  necessary  to  obtain  the  missing  values  from  other  investigators. 

It  should  be  noted  that  the  SEMS  acquired  surface  tension  data  at 
approximately  three-minute  intervals.  Only  periods  with  generally 
acceptable  surface  tension  values  are  listed  in  the  catalog. 

About  5 /  of  these  data  are  suspicious  and  require  additional 
examination  because  of  possible  artifacts  associated  with  the 
measurement  system  operation. 
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EG  INNING  10-20-1503 


Figure  7.  Natural  Slicks,  Surface  Tension  - 


After  final  processing,  the  surface  tension  data  will  be 
corelated  with  radar  backscatter  results  and  with  internal  wave 
events  acquired  simultaneously  during  the  TOWARD  Experiment. 
These  data  intercomparisons  are  of  interest,  because  internal 
waves  play  a  role  in  the  formation  of  coherent  banded  organic 
surface  films,  while  the  organic  films  in  turn  attenuate  and 
resist  the  formation  of  capillary  waves  and  reduce  microwave 
radar  backscatter. 
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APPENDIX  I 


COMPUTER  PROGRAM  FLOWCHART 


I  •  Surface  Tension 


MOVE  MEMORY 
CALL  PART  II 


The  location  of  BASIC  programs  is  changed 
to  reserve  memory  for  the  graphics.  Part 
II  is  called. 


Surface  Tension  II 


An  initialization  step. 


Pump  begins  flushing  system  early  during 
the  first  on  cycle. 


When  the  system  is  first  turned  on  the 
date  and  time  must  be  entered. 


Printer  labels  paper.  Graph  axes  are 
drawn  on  screen.  Disk  file  is  automatically 
named  with  date  and  time.  (Repeated  hourly). 


Time  and  sensor  signal  are  shown  on 
screen  while  water  sample  fills  sample 
chamber  and  overflows. 


The  number  of  points  determines  how 
long  water  is  pumped  through  the  sample 
chamber. 


The  pump  shut-off  time  is  recorded  as  the 
sample  time. 


As  water  drains  slowly  from  the  chamber 
the  signal  is  read  rapidly  400  times. 
Values  are  stored. 


The  maximum  value  occurs  when  the  glass 
plate  detaches  from  the  water  surface. 


Points  prior  to  the  maximum  and  10  points 
after  it  are  ignored.  The  remaining 
points  of  the  400  are  averaged  for  the 
minimum. 
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COMPUTE 

SURFACE 

TENSION 


Surface  tension  =  (max  -  min)  X 
(calibration  factor). 


r 


PRINT  TIME 
&  SURFACE 
TENSION 


Date,  time,  surface  tension,  max,  min, 
and  difference  are  printed  on  the  hard 
copy . 


PLOT  POINT 
ON  SCREEN 


I  The  cycle  continues  until 
|  stopped  with  the  ESC  key. 


Points  are  taken  approximately  every 
3  minutes. 


At  the  end  of  an  hour  the  data  is  saved 
on  disk,  the  screen  graph  is  drawn  on 
the  hard  copy,  a  new  file  is  created, 
the  screen  is  erased  and  new  axes  are 
drawn,  etc. 


SAVE  OATA 
ON  OISK 
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PROJECT  1UUAHD 

URL  Surface  Tension  Measurement  System 

I 


Catalog  of  Surface  Tension  Data 


Date 

T  ime 

Date 

T  ime 

09-20-84 

1620-1648 

10-15-84 

1616-1703 

09-21-84 

1624-1657 

10-16-84 

0933-1 146 
1157-1759 

09-22-84 

0836-0851 

1833-1842 

0915-1013 

10-17-84 

0852-1414 

09-24-84 

0822-0958 

101 3-1057 

10-18-84 

0915-1238 

1105-1232 

1315-1357 

1240-1358 

1416-1440 

09-25-84 

0814-0937 

10-20-84 

0940-1144 

1014-1058 

1315-1815 

1109-1333 

1853-2223 

09-26-84 

0853-0856 

10-23-84 

0843-1006 

0907-1459 

1  152-1422 
1634-1722 

09-27-84 

0847-1459 

10-24-84 

0847-1128 

09-28-84 

0839-1210 

1208-1448 

1223-1501 

1516-1654 

10-02-84 

0924-1044 

10-25-84 

0835-1256 

1053-1558 

1454-1557 

10-03-84 

0906-1520 

10-27-84 

0847-1021 

1  124-1654 

10-04-84 

0842-0957 

10-28-84 

1110  ’250 

10-06-84 

1005-1150 

1129-1609 

1315-1558 

10-29-84 

1109-1619 

10-09-84 

0556-0659 
0850-091  8 

10-30-84 

0823-1056 

1341-1359 

1424-1450 

1120-1519 

10-31-84 

0936-1256 

10-11-84 

1843-2008 

2030-2105 

1415-1457 

10-12-84 

0841-1317 

1651-1814 
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WAVE  FOLLOWER  MEASUREMENTS  DURING  TOWARD  84/85 

Paul  A.  Hwang 

Ocean  Research  and  Engineering 
La  Canada,  California  91011 


FOREWORD 


As  a  part  of  the  TOWARD  experiment,  the  wave  follower 
measurements  were  to  provide,  among  other  things,  ground  truth 
of  short  waves  for  radar  backscatter  and  SAR  imaging  groups. 

This  report  presents  a  short  account  on  the  scope  of  the  measure¬ 
ments  and  the  progress  of  the  data  analysis.  The  report  is 
intended  for  internal  use  of  the  TOWARD  investigators.  The 
purpose  is  to  specify  the  available  data  sets,  define  the 
analysis  in  progress  and  provide  up-to-date  results  to  potential 
TOWARD  users. 
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1.  INSTALLATION  AND  INSTRUMENTATION 


The  wave  follower  assembly  was  installed  on  the  south  extension  structure 
of  the  NOSC  Tower  (see  Figure  1).  The  assembly  was  placed  20  feet  away  from 
the  tower  structure  and  connected  to  it  by  a  catwalk.  The  details  of  the  wave 
follower  and  instrumentation  are  presented  elsewhere  (Shemdin,  1980;  Shemdin 
and  Hoff,  1986).  The  following  gives  a  brief  description  of  the  system. 

A  set  of  instruments  were  installed  at  the  bottom  of  the  fixed  structure 
(see  Figure  2)  to  provide  long  wave  directional  measurements.  For  this  pur¬ 
pose  a  pressure  transducer  (a)  and  a  two-dimensional  E-M  current  meter  (b)  were 
used. 

On  the  moving  frame  (also  referred  to  as  C-frame),  four  instruments  were 
mounted,  as  shown  in  Figure  3.  Of  these  a  capacitance  gage  was  used  as  an 
elevation  control  sensor  (c).  This  sensor  provided  the  feedback  signal  to  the 
servo-system  to  maintain  the  C-frame  at  a  fixed  elevation  with  respect  to  the 
instantaneous  water  level.  A  two-dimensional  E-M  current  meter  (d)  was  mounted 
on  the  lower  arm  of  the  C-frame  to  provide  the  two  horizontal  components  of 
the  current.  The  two-dimensional  optical  wave  slope  sensor  (e)  and  the  optical 
wave  height  sensor  (Ret icon  camera,  (f)  were  both  used  for  measuring  high  fre¬ 
quency  short  waves.  Both  optical  sensors  used  the  10  mW  laser  (g)  that  was 
mounted  on  the  lower  arm  of  the  C-frame  for  a  light  source. 

The  wave  slope  sensor  was  calibrated  in  the  laboratory  before  and  after 
each  experiment.  Light  refraction  at  the  air-glass  interface  was  used  to 
simulate  refraction  of  light  at  the  air-water  interface  (see  top  insert  of 
Figure  4).  The  magnitude  of  the  refracted  angle,  was  varied  by  moving  the 
glass  lens  horizontally.  By  rotating  the  lens  on  the  horizontal  plane,  changes 
in  azimuthal  angle,  ^  were  achieved  (see  bottom  insert  in  Figure  4).  A 
typical  calibration  result  is  shown  in  Figure  4.  Averaged  over  the  azimuthal 
angles.  Figures  5a,  b  show  the  calibration  curves  for  phases  I  and  II  of  the 
experiment,  respectively.  The  uncertainty  is  estimated  to  be  the  lesser  of 
±0.5°  or  ±13%  in  the  range  0  to  30°  in  water  surface  slope.  This  result 
is  derived  from  phase  II  calibration  tests  where  eight  repetitions  were 
performed.  Only  two  repetitions  were  performed  in  phase  I  calibration  tests. 
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Figure  2.  Wave  Follower  Structural  Assembly, 


Figure  3.  Wave  Follower  Instrument  Frame.  Capacitance  Gage 
is  used  as  Elevation  Control  Sensor  (c) ;  Two-dimensional 
current  meter  is  used  for  measuring  two  horizontal  currents 
components  below  water  surface  (d);  Optical  Receiver  records 
two  slope  components  of  surface  normal  (e) ;  Reticon  camera 
records  vertical  displacement  of  water  surface  (f)  and  using 
He-Ne  Laser  ( g)  . 
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Laboratory  Calibration  of  Wave  Slope  Instrument  Used  in  TOWARD  84/85. 
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2.  DATA  REDUCTION 


During  the  two  phases  of  the  TOWARD  experiment,  more  than  130  hours  of 
data  were  recorded,  from  which  20  hours  of  data  were  selected  for  further 
studies.  The  criteria  of  selection  were  based  on  periods  of  aircraft  flights, 
periods  of  stereo-photo  sessions  and  range  of  wind  and  wave  conditions.  The 
environmental  conditions  varied  from  steady  (for  periods  of  at  least  30 
minutes)  to  transient  (either  in  wind  speed  or  wind  direction).  The  waves 
ranged  from  large  swell  to  strong  local  sea.  Tdbles  I  and  II  list  the 
conditions  of  the  data  segments  that  were  digitized.  At  the  time  of  this  data 
compilation,  the  swell  information  had  not  been  completed;  it  will  be  included 
in  the  final  report.  In  Figure  6  the  distribution  of  data  segments  versus 
wind  speed  is  shown.  In  phase  II,  two  storms  swept  through  the  tower  region 
which  generated  wind  speeds  up  to  23  knots  (12.5  m/s).  The  latter  events 
broadened  significantly  the  range  of  wind  speed  measured. 


Table  I.  Data  Segments  Recorded  in  Phase  I. 


Date 

Pacific* 

Time 

Wind 

Speed  (m/s) 

Wind 

Direction 

Stereo 

Photo 

Flights 

10/17/84 

1315  -  1345 

7.5 

285 

CV990 

10/20/84 

1355  -  1455 

6 

265 

2355  -  0025 

4 

120 

10/29/84 

1435  -  1505 

5 

300 

1520  -  1550 

4.5 

300 

10/30/84 

1535  -  1605 

3 

290 

10/31/84 

1030  -  1100 

2.5 

235 

1300  -  1600 

3-5 

280-300 

4  Series 

CV990 

F4 

11/03/84 

1820  -  2045 

2.5 

270 

2330  -  2400 

3 

25 

11/04/84 

0030  -  0130 

3 

10 

0945  -  1230 

2 

220 

1  Series 

F4 

1530  -  1630 

3 

300 

2  Series 

1835  -  1930 

1 

175 

2320  -  2350 

2.5 

0 

11/07/84 

1625  -  1655 

4 

260 

3  Series 

CV990 

*Pac if ic 

Time  is  8  hours 

lag  Greenwich 

Time . 

Table  11.  Data  Segments  Recorded  in  Phase  II. 


Date 

Time 

Wind 

Speed  (m/s) 

Wind 

Direction  (°) 

Stereo 

Photo 

Flights 

3/19/85 

1400  - 

1530 

6 

260 

CV990 

F4,  OV-ID 

3/20/85 

1000  - 

1130 

6 

330 

3/21/85 

1450 

•  1500 

3 

310 

1  Series 

F4 ,  OV-ID 

2030  - 

2230 

0 

/ 

3/22/85 

0340  - 

-  0410 

1.5 

90 

1145  - 

1215 

3 

290 

2  Series 

3/24/85 

1330  - 

•  1400 

3 

270 

3/25/85 

1120  - 

■  1150 

5 

235 

OV-ID 

1340 

-  1700 

7 

240 

1730  - 

1800 

6 

230 

1940 

-  2010 

2.5 

180 

3/26/85 

1011  - 

-  1140 

4 

210 

4  Series 

F4 ,  OV-ID 

3/27/85 

1020  - 

-  1050 

8 

225 

1  Series 

1330  - 

-  1530 

9± 

225 

CV990 

3/28/85 

0930 

-  1130 

12± 

225 

2  Series 

4/1/85 

2100  - 

-  2230 

0 

3.  WAVE  SLOPE  DATA  ANALYSIS 


The  SAR  committee  assigned  "priority"  to  the  L-band  SAR  flight  of 

October  31,  1984,  1300  to  1500  hours  local  time.  The  wave  follower  effort  was 

directed  to  study  this  particular  data  set.  A  sample  of  the  capillary  wave 

slope  data  are  shown  in  Figure  7.  Here,  nx>  denotes  slope  in  the 

upwind-downwind  direction,  and  denotes  slope  in  the  cross-wind 

direction.  The  surface  displacement  of  long  waves,  n»  is  measured  by  a 

submerged  pressure  gage  and  is  also  shown  in  Figure  7.  The  starting  time  for 

this  data  series  is  1300  local  time.  The  wind  speed  at  this  time  was  3.5  m/ s 

and  the  wind  direction  was  280°.  The  swell  period  was  8  to  10  seconds,  from 

direction  290°.  The  swell  significant  wave  height  was  1.2m.  The  variance  of 

2  -3 

slopes  in  the  two  directions  was  found  to  be  very  similar  ( n  *  3.72  x  10 
2-3 

and  n  *  3.78  x  10  ).  In  the  following  the  results  of  analysis,  up 

to  date,  are  described. 
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Figure  6.  Experimental  Conditions  in  TOWARD  84/85. 
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ire  7.  Records  of  Capillary  Wave  Slopes  in  the  Upwind-Downwind 
Cross-Wind  Directions,  and  Water  Surface  Displacement  Measured 
.  Submerged  Pressure  Gage. 


3.1  The  Probability  Density  Distribution  of  Wave  Slopes 


The  probability  density  functions,  PDF's,  of  the  two  slope  signals  are 
found  to  be  somewhat  different,  as  shown  in  Figure  8a.  As  observed  in 
previous  studies  of  sun  glitter  (Cox  and  Monk,  1954),  the  upwind-downwind 
component  of  slope  is  more  skewed  compared  to  the  cross-wind  components.  For 
this  wave  record,  the  skewnesses,  defined  as 


and 


(1) 


(2) 


were  found  to  be  -2.04  and  -0.60  for  the  x-  and  y-components,  respectively. 

In  Equations  (1)  and  (2)  the  means  of  the  signals  are  zero  and  the  upper  bar 
indicates  time  average.  At  low  wind  speeds,  such  as  prevailed  in  this  record, 
the  PDF's  are  typically  narrower  than  normal  distribution.  The  peakness  of 
distribution  is  defined  as 


and 


(3) 


(4) 


They  were  calculated  for  these  records  to  be  10.94  and  19.02,  respectively. 
The  normalization  factors  in  Equations  (1)  through  (4)  are  related  to  the 
variances  of  the  two  slope  components,  which  as  indicated  before,  are  very 
close  in  value.  The  standard  deviation  is  equal  to  the  rms  slope  which 
becomes  0.061  (or  approximately  3.5°)  for  slopes  in  both  directions.  From 
Figure  8,  it  is  obvious  that  most  waves  have  slopes  less  than  or  equal  to 
10°.  There  are,  however,  long  tails  on  both  end  of  the  PDF  curves.  Here 
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slopes  as  high  as  28°  (excluding  "dropouts"  from  the  measurements  -  discussed 
later  in  the  report)  were  registered.  There  is  considerable  interest  in  these 
high  slope  tails  as  they  appear  to  correlate  with  the  specular  component  of 
radar  backscatter  signals. 

In  the  subsequent  analysis,  we  apply  frequency  filters  to  the  slope 
signals.  Figures  8b,  c,  and  d  show  the  PDF's  for  the  2.5  Hz  low-pass  filter, 
2*5  -  7.5  Hz  band  pass  filter  and  7.5  Hz  high  pass  filter,  respectively.  For 
comparison,  the  rms  slopes  used  for  normalization  of  Figures  8b,  c  and  d  were 
those  used  for  the  unfiltered  signals  (Figure  8a).  Comparison  of  these  three 
plots  indicates  that  the  skewnesses  of  PDF's  are  mainly  contributed  by  the  low 
frequency  range,  or  frequencies  less  than  2.5  Hz  in  the  wave  spectra.  The 
steepnesses  of  these  waves  are  generally  small  (less  than  7°).  Waves  in  the 
middle  frequency  range  (2.5  -  7.5  Hz)  show  more  symmetric  slope  distributions 
(see  Figure  8c).  Here  the  wave  slopes  are  still  low  (less  than  7°).  The 
similarity  of  the  upwind-downwind  and  cross-wind  wave  slopes  in  this  middle 
frequency  range  indicates  a  more  isotropic  distribution  of  short  waves. 
Finally,  from  Figure  8d,  it  is  clearly  shown  that  the  steep  components  of  wave 
slopes  come  entirely  from  the  high  frequency  waves.  The  characteristics  of 
long  tails  and  narrow  distributions  in  Figure  8d  are  similar  to  those  of  the 
unfiltered  signals,  Figure  8a. 

The  PDF's  in  the  mid-  and  low-frequency  ranges  are  very  close  to  Gaussian 
when  normalized  by  their  respective  rms  slopes,  as  shown  in  Figures  9b  and  c. 
The  distribution  of  the  high  frequency  waves  is  very  narrow  (see  Figure  9d). 
The  unfiltered  PDF  is  shown  in  Figure  9a  for  comparison.  Table  III  shows  the 
partitioning  of  variance  in  these  three  frequency  ranges  as  well  as  the 
skewness  and  flatness  (kurtosis)  statistics.  From  this  table  the  following 
observations  are  noted: 

2  2 

(a)  The  ratio  of  q  /q  varies,  2.17  at  low  frequency,  1.29  at  mid-range, 

x  y 

and  0.65  at  high  frequency.  This  reflects  the  gradual  change  from 
two-dimensional  (organized)  long  waves  to  a  more  random  short  wave 
distribution.  The  ratio  0.65,  at  the  high  frequency  end  is  somewhat 
puzzling.  The  time  series  in  Figure  7,  (also,  the  drop-out 
statistics  which  will  be  discussed  later)  suggest  that  the 
cross-wind  component  of  wave  slope  may  contain  more  noise  compared 
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Figure  8a.  PDF  of  Wave  Slopes,  Unfiltered.  The  RMS  Variance  of 
Unfiltered  Signal  Is  Used  for  Normalization  of  Standard  Deviation 
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Figure  8c.  PDF  of  Wave  Slopes,  2.5  to  7.5  Hz  Band  Passed 
Filtered.  The  RMS  Variance  of  Unfiltered  Signal  Is  Used 
for  Normalization  of  Standard  Deviation. 
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Figure  8d .  PDF  of  Wave  Slopes,  7.5  Hz  High  Pass  Filtered  Results.  The 
RMS  Variance  of  Unfiltered  Signal  Is  Used  for  Normalization  of  Standard 
Deviation . 
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Figure  9a.  PDF  of  Wave  Slopes,  Normalized  by  the 
Total  Variance,  Unfiltered. 
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Figure  9c.  PDF  of  Wave  Slopes,  Normalized  by  the  Filtered  Variance, 
2.5  to  7.5  Hz  Band  Pass  Filtered. 
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Figure  9d .  PDF  of  Wave  Slopes,  Normalized  by  the  Filtered  Variance, 
7.5  Hz  High  Pass  Filtered  Results. 
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to  the  upwind-downwind  channel.  A  more  realistic  ratio  would  be  1 
for  an  isotropic  short  wave  field.  More  analysis  of  data  is  needed 
to  confirm  this  observation,  and  is  planned. 
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(b)  The  skewness  of  the  filtered  signal  is  generally  smaller  than  the 
unfiltered  signal  except  for  the  high  frequency  component  in  the 
cross-wind  direction.  The  small  magnitudes  of  skewness  in  the  low 
to  intermediate  frequency  bands  of  the  cross-wind  slopes  indicate 
possible  two-dimensionality  of  these -waves ,  i.e.,  long-crested  waves. 

(c)  The  flatness  (kurtosis)  estimates  show  an  increasing  trend  toward 
the  high  frequency  wave  components.  At  the  low  frequency  end,  the 
magnitude  is  close  to  that  of  a  Gaussian  distribution  (for  which, 
the  kurtosis  is  3).  An  increase  above  this  value  indicates  a 
tendency  toward  a  narrower-band  frequency  distribution.  The  larger 
the  magnitude  of  this  statistical  measure,  the  more  the  occurrences 
of  very  large  and  very  small  wave  slopes,  compared  to  a  Gaussian 
distribution. 

Cox  and  Munk  (1954)  found  from  their  sun  glitter  measurements  that  the 
skewness  in  the  cross-wind  direction  (their  =  0.11  ±  0.03)  is  less 
than  the  upwind-downwind  (Cq-j  =  -0.42  ±  0.12)  at  high  wind  speed  (14  m/s); 
and  relatively  small  (nearly  zero)  at  very  low  winds  (velocity  unspecified). 
The  flatness  (peakness  in  their  terminology)  estimates  from  their  measurements 
shows  higher  cross-wind  (C^  *  0.4  ±  0.2)  than  upwind-downwind  (C^  = 

0.2  ±  0.4)  components.  These  results  are  consistent  with  our  present 
measurements.  Our  frequency  decomposition  analysis,  however,  indicates  that 
the  distribution  of  wave  slopes  at  different  frequency  bands  are  quite 
different  from  each  other.  More  detailed  studies  on  wind  speed  and  swell 
dependence,  and  of  longer  time  series,  are  forthcoming. 

3 . 2  Longitudinal  Variation  of  the  Slope  Variances 

For  the  SAR  and  radar  return  studies,  the  longitudinal  (or  spatial) 

variation  of  short  waves  along  the  long  wave  profile  is  of  considerable 

interest.  The  period  of  the  swell  on  October  31,_1984  v?as  in_the  range  8  to 

2  2  2 

10  seconds.  In  Figure  10  the  mean  square  slope  Vq  ”  ( n  ♦  q  ),  calculated 

x  y 


tv 
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e  10.  Longitudinal  Variation  of  Mean  Square  Slope  of  Short  Waves 
ata  Dropouts  Along  Profile  of  Long  Waves. 


from  one  second  record  length,  is  plotted  in  relation  to  the  long  wave 
profile.  The  large  variation  of  the  short  waves  at  different  phases  of  the 
long  wave  is  clearly  demonstrated.  The  maximum  and  minimum  of  the  mean  square 

-2  -3 

slopes  (one-second-averages)  are  found  to  be  1.42  x  10  and  2.23  x  10  , 

respectively.  The  mean  value  is  7.5  x  10  .  The  fluctuation  is  nearly  100% 

of  the  mean.  When  the  small  steepness  of  the  long  wave  is  considered  (Hg  = 
1.23m,  Tg  =  9s,  the  "average"  slope  is  approximately  0.01)  the  modulation  of 
short  waves  is  approximately  100  times  that  of  long  wave  slope.  Such  a  high 
degree  of  modulation  is  indeed  significant.  The  local  maxima  of  the  slope 
variance  generally  occur  near  wave  crests.  It  is  noted  that  secondary  maxima 
occur  (not  necessarily  of  smaller  magnitudes)  between  the  long  wave  crest. 

The  peaks  that  appear  above  the  wave  troughs  are  located  at  12  and  20  seconds 
along  the  time  record.  The  strong  local  disturbances  and  secondary  peaks 
suggest  nonlinearity  of  short  waves  and  possibly  instability  along  the  long 
waves.  This  observation  of  large  modulation  of  short  waves  by  the 
low-amplitude  long  waves  is  difficult  to  explain  by  a  two-scale  model  such  as 
that  proposed  by  Longuet-Higgins  and  Stewart  (1960)  and  Longuet-Higgins  (1969). 

In  Figure  10  the  number  of  data  points  that  have  a  slope  magnitude  0.5  or 
greater  (this  value  is  close  to  the  theoretical  wave  breaking  limit)  are  shown 
in  relation  to  the  long  wave  profile.  These  points  are  denoted  as 
"dropouts."  They  bear  a  relationship  to  breaking  waves.  Note  that  at  the 
time  of  this  writing,  the  dropout  count  includes  those  points  that  have  a 
negative  slope  less  than  -0.5.  For  correlation  with  breaking  waves  these 
points  should  be  excluded.  The  mean  square  slope  results  shown  here  do  not 
include  the  dropout  points.  The  raw  data  had  been  edited  by  replacing  the 
dropout  values  by  the  mean  value  of  those  points  before  and  after  the  dropout 
group.  The  largest  single  group  of  dropouts  occurs  between  14  and  15 
seconds.  Here,  a  consecutive  fourteen  data  points  exceeded  the  threshold 
magnitude . 

3 . 3  Wave  Slope  Spectra 

3.3.1  Measured  Frequency  Spectra 

The  frequency  spectra,  shown  in  Figure  11,  give  the  energy  distribution 
of  wave  slopes  at  various  frequencies.  Both  the  upwind-downwind  and 
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cross-wind  slope  spectra  are  included.  The  slopes  of  the  spectra  at  the  high 
frequency  ends  are  both  approximately  -1.2.  Similar  to  what  was  observed  from 
the  probability  distribution  of  waves  (Section  3.1)  at  low  frequency,  the 
spectral  energy  level  in  the  upwind-downwind  direction  is  much  higher  than  in 
the  cross-wind  direction.  The  two  curves  gradually  merge  together  at  a 
frequency  close  to  10  Hz.  There  is  a  sudden  step  up  of  energy  level  in  the 
cross-wind  component  at  23  Hz.  A  plateau  is  formed  which  extends  to  43  Hz. 
Beyond  this  frequency  the  spectrums  follows  a  -1.2  slope  again.  The  region  of 
the  plateau  and  the  elevated  energy  level  in  the  cross-wind  signal  cannot  be 
properly  explained  at  the  present  time. 

The  frequency  spectra  computed  from  one-second  data  segments  are  shown  in 
Figure  12  in  relation  to  phase  location  along  the  long  waves.  Different 
shapes  are  observed  at  various  phase  locations  along  the  long  wave  profile. 
Similar  to  the  distribution  of  mean  square  slopes,  such  variation  demonstrates 
the  modulation  of  short  waves  by  the  long  waves.  Inspection  of  each 
individual  spectrum  indicates  that  near  the  trough  regions,  the 
upwind-downwind  spectral  slopes  are  generally  smaller  (flatter),  with  slopes 
ranging  from  -0.6  to  -0.2,  compared  to  other  regions.  The  combination  of  the 
low  mean  square  slope  and  a  broad  spectrum  at  the  trough  region  implies  strong 
interaction  between  wave-induced  current  and  the  short  waves. 

3.3.2  Conversion  to  Wave  Number  Spectra 

For  the  SAR  imaging  and  radar  backscatter  applications,  wave  number 
(versus  wave  frequency)  spectra  are  more  useful.  Conversion  of  wave 
information  either  from  time  domain  to  space  domain  (Evans  and  Shemdin,  1980), 
or  from  frequency  domain  to  wave  number  domain  (Hughes,  1978;  Atakturk,  1984) 
is  not  a  simple  task.  After  reviewing  the  various  methods  reported,  the 
method  employed  by  Hughes  is  adopted  for  this  study.  The  method  was  discussed 
in  detail  in  his  paper  of  1978.  Briefly,  the  spectral  energy  at  a  particular 
frequency  G( o)  da  has  its  contributions  from  all  wave  numbers  and  wave 
directions  that  satisfy  the  Doppler  shifted  dispersion  relationship. 


TIME  (s> 

Figure! 2.  Modulation  of  Frequency  Spectra  Along  The  Long  Wave, 
(a)  Slope  Spectra  of  Cross-Wind  Component,  (b)  Slope  Spectra 
of  Upwind-Downwind  Component,  (c)  The  Long  Wave  Profile. 


and  can  be  written  as 


G(o)  do  = 


// 


F(k)D( e)  B(k,9)  5(0,  k,  u,  e)  dedk  , 


(6) 


where  o  is  the  encounter  frequency  of  the  short  waves,  og  is  the  intrinsic 
frequency,  explicitly  expressed  as 

2  _  .  .  .3  , 

o  -  gk  +  tk  ,  • 

k  is  the  wave  number  vector  with  magnitude  k  and  direction  9,  u  is  the 
current  vector  with  magnitude  u  and  direction  a,  F(k)D(g)  is  the  wave 
number  spectrum  with  k  and  0  dependence  assumed  separable,  B(k,9)  the 
weighting  function,  its  form  depends  on  the  type  of  signal  analyzed.  More 
explicitly 

1  for  wave  height  spectrum  (S  ) 


(7) 


2  2 

k  cos  e  for  upwwind/downwind  slope  spectrum  (S 
2  2 

B(k,e)  *  k  sin  9  for  cross-wind  slope  spectrum  (S  ) 

V 

k  cos  0  for  height,  upwind-downwind  and  cross-spectrum  (S  ) 


nn. 


k  sin  9  for  height  and  cross-wind  cross-spectrum  (S  ), 


nn 


6(0,  k,  u,  0)  is  the  delta  function  such  that  when  k  satisfies  the 
dispersion  relation  (Equation  5),  the  value  is  1,  otherwise  it  is  0;  g  is  the 
gravitational  acceleration;  and  t  the  surface  tension  of  water. 


Equation  6  can  be  written  in  matrix  form 

g  =  A  f  , 

where  £  is  the  vector  of  which  the  elements  are  frequency  spectral  density 
G(o),  £  the  corresponding  wave  number  spectral  density  F(k^),  matrix  A 


(8) 


connects  £  and  £  and  with  elements  expressed  as 

a..  =  7  D( 9 . )  B(k.,0.)  k.  A9. 

ij  Lu  J  J  J  J  J 


(k , 9)  • 


3k 

3o 


(9) 


ij 
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where  the  (k,©)^  pairs  are  those  that  satisfy  dispersion  relation  for  a 
particular  frequency  o,  and  current  u.  Note  that  in  Equation  9,  the 
absolute  value  of  3k/3o  was  used  since  do  on  the  left-hand  side  of 
Equation  6  is  the  frequency  "width"  such  that  G( a)  do  is  the  spectral 
energy  at  frequency  o,  and  kdSdk  is  the  integration  "area"  such  that  the 
expression  on  the  right-hand  side  of  Equation  6  denotes  partitioning  of 
G(a)  do  in  different  wave  numbers.  Both  the  "width"  and  "area"  need  to  be 
positive. 

The  algorithm  developed  was  tested  with  a  hypothetical  spectrum  and  found 
to  be  satisfactory  (details  in  the  Appendix).  In  application  to  the  measured 
data,  the  frequency  spectra  calculated  from  one-second  records  (Figure  12) 
were  substituted  for  £  vector  in  Equation  8,  the  matrix  A  was  readily 
constructed  with  a  given  directional  distribution  of  the  form 


D( e)  =  cos  0 


<  a  <  e2 


(10) 


In  the  following  presentation,  p  =  2,  0^  =  -ir/2,  0^  =  ir/2.  The  necessity  of 

prescribing  D( @)  is  believed  not  to  be  a  serious  handicap  since  the 

measurements  of  nx  and  ry  provides  partial  information  of  wave 

directional  distribution.  Possible  forms  of  D(e)  were  suggested  by  Hughes 

(1978)  and  Phillips  (1985)  among  others.  In  the  TOWARD  Experiment, 

stereo-photo  pairs  were  taken  at  different  phase  locations  along  the  long 

wave.  The  stereo-photos  are  currently  being  analyzed.  The  results  will 

clarify  much  of  the  uncertainties  on  the  directional  distribution  of  short 

waves.  The  inversion  scheme  developed  is  flexible  and  can  accommodate  any 

form  of  D( 0) .  For  a  preliminary  computation  Figure  13  presents  the  results 

of  spectral  conversion  from  frequency  to  wave  number.  The  wave  number  spectra 

(of  wave  height)  F(k)  showed  more  similarity  in  spectral  shape  along  the  long 

wave  profile  compared  to  its  frequency  counterpart.  The  39  consecutive  wave 

number  spectra  (only  21  presented  in  Figure  13)  all  can  be  fitted  to  form 
-4 

k  ,  with  the  exception  of  1,  2,  10,  11,  19,  and  20.  The  wave  number 

-4 

spectra  that  deviate  from  k  occur  exclusively  in  the  trough  regions  of 
long  waves  and  appear  to  have  low  mean  square  slopes.  For  some  wave  troughs 
(segments  10,  11  and  19,  20)  the  wave  number  spectra  still  have  the  k  ^ 
shape  at  the  short  wave  range,  but  the  energy  contents  of  the  longer  waves  are 
considerably  less.  The  fact  that  nearly  all  the  wave  number  spectra  are 
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saturated  or  super-saturated  seems  to  imply  that  short  waves  have  small 
relaxation  times.  This  is  in  agreement  with  the  laboratory  studies  of  Wu 
(1979)  on  decay  and  regeneration  of  ripples. 

The  wave  number  spectra  converted  from  the  cross-wind  slope  are  less 

uniform,  but  still  show  consistent  agreements  with  those  obtained  from  the 

upwind-downwind  slopes.  The  averaged  wave  number  spectrum  over  the  entire 
•  "“A 

profile  of  long  waves  also  shows  the  k  dependence,  as  shown  in  Figure  14. 

! 

4.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

In  this  report,  we  have  presented  the  analysis  of  short  wave  modulation 
by  long  waves  based  on  the  two-dimensional  wave  slope  measurements  conducted 
in  the  TOWARD  Experiment.  At  the  time  the  present  record  was  acquired 
(10-31-84,  1300  Pacific  Time),  the  wind  speed  was  3.5  m/s,  from  direction 
280°.  The  swell  recorded  had  a  period  in  the  range  8  to  10  seconds  and  a 
significant  wave  height  of  1.2  m.  The  swell  direction  was  from  direction 
290°.  During  this  period  the  wave  follower  performed  exceptionally  well. 
Based  on  this  criterion  of  data  dropout  (tane  £  0.5),  the  data  quality  in 
the  upwind-downwind  direction  is  considered  excellent  (with  dropout  rate  of 
0.47%.  The  cross-wind  slope  data  quality  is  considered  good  (3.49%  dropout 
rate).  The  accuracy  of  the  slope  measurements  was  ±0.5°  at  small  slopes  and 
±13%  of  steepness  near  30°  (Figure  4). 

The  results  from  this  analysis  shows  that: 

(a)  The  probability  distribution  of  the  upwind-downwind  slope  is 

skewed.  Tne  most  probable  surface  slope  is  approximately 

0.7  degrees.  Although  the  distribution  is  bell  shaped,  is  is  much 

narrower  than  Gaussian  with  a  flatness  factor  of  10.94.  The  PDF  of 

cross-wind  slope  is  more  symmetric  and  narrower  (see  Table  III). 

When  frequency  filters  are  applied,  it  is  obvious  that  the  skewness 

comes  from  the  low  frequency  components  while  peakness  comes  from 

the  high  frequency  components.  Furthermore,  the  large  steepness  is 

derived  exclusively  from  waves  with  frequencies  higher  than  7.5  Hz. 

S lope s  of  waves  below  7.5  Hz  seldom  exceed  7°.  The  ratio  of 

2  2.. 

n  /n  for  the  unfiltered  signals  is  close  to  1,  but  differs 
x  y 

considerably  in  the  different  frequency  regions.  The  range  is  from 
2.17  at  low  frequencies  (£2.5  Hz),  1.29  at  mid-range  (2.5  -  7.5  Hz) 


. I 


0 


1 


LOG  [k  (1/cm)] 

Figure  14,  The  Phase  Averaged  Wave  Number  Spectrum  of  39  Segments 
(Each  of  One  Second  Data  Length)  for  the  Upwind-Downwind  Component 
of  Wave  Slope  Measurement. 


and  0,65  at  high  frequency  range  (<  7.5  Hz).  The  PDF  distribution 
is  closest  to  Gaussian  at  low  frequency,  and  most  symmetric  in  mid 
range.  The  similarity  of  the  normalized  PDF  of  the  two  slope 
components  at  raid-range  (wave  lengths  from  X-band  to  L-band)  is 
illuminating  (Figure  8c). 

(b)  The  frequency  spectra  averaged  along  the  long  wave  profile  display 

-1.2. 

the  form  f  ‘  in  both  upwind-downwind  and  cross-wind  directions 

(Figure  10).  This  implies  a  close  to  saturation  condition  even  at 

3.5  m/ s  wind  speed  condition.  Similarly,  the  averaged  wave  number 

-4 

spectra  show  the  form  k  (Figure  13). 

(c)  Strong  modulation  of  short  waves  by  long  waves  (Figure  9).  The 
magnitude  of  modulation  is  nearly  100  times  of  the  long  wave  slope. 

(d)  The  frequency  spectra  are  not  only  subject  to  the  amplitude 
modulation  similar  to  the  mean  square  slopes,  but  also  show 
significant  variation  in  the  spectral  slopes  or  equivalently  the 
spectral  band  width.  The  spectra  are  considerably  broader  over  the 
trough  regions  where  the  wave  induced  orbital  velocity  is  in  the 
direction  opposite  to  wind  direction.  The  broadening  of  spectra  in 
an  adverse  current  region  implies  active  interaction  between  the 
wave  induced  current  and  the  short  waves. 

(e)  The  observation  in  (d)  are  also  seen  in  the  converted  wave  number 

spectra  along  the  long  wave  profile.  While  at  most  phase  locations 

-4  ... 

along  the  long  waves  profile  a  k  relationship  is  observed, 
deviations  from  such  a  form  is  observed  in  the  trough  regions.  The 
appearance  of  a  saturated  or  supersaturated  spectrum  suggests  that 
fast  regeneration  of  modulated  short  waves  is  effected  by  wind. 

The  results  discussed  are  based  on  a  relatively  short  record 
(39  seconds).  Subsequent  computations  on  longer  time  series  did  not  change 
the  above  observations,  except  in  minor  details.  In  the  continuing  data 
analysis,  we  expect  to  develop  schemes  for  ensemble  averaging  for  better 
statistical  representation,  and  also  to  include  various  environmental 
conditions . 

There  is  no  doubt  that  more  information  can  be  extracted  from  the 
wave-follower  data  sets  acquired  in  the  TOWARD  Experiment. 
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APPENDIX  -  CONVERSION  OF  FREQUENCY  TO  WAVE  NUMBER  SPECTRA 


The  scheme  for  the  spectral  conversion  from  frequency  to  wave  number 

-4 

domain  described  in  the  paper  was  tested  using  a  hypothetical  k  spectrum. 

The  frequency  spectra  for  different  current  velocities  were  calculated 

.  2 
directly  from  Equation  6  with  D ( 0 )  =  cos  0.  Two  ranges  of  0,  ±  ir/4  and 

±  ir/2,  were  used,  but  no  significant  difference  was  found.  The  results 

presented  in  the  following  is  for  -  ir/4  <  0  <  tt/4. 

The  treatment  of  the  dispersion  relationship  (Equations  5  and  7)  is 
slightly  different  from  that  of  Hughes  (1979).  In  the  original  design, 

Hughes,  based  on  the  symmetricity  of  spectrum,  chose  to  replace  on  the  left 
hand  side  with  ±  a,  and  adopted  both  roots  of  the  right  hand  expression  in 
Equation  7  (here,  surface  tension  is  set  to  0).  This  approach  makes  no 
difference  if  u  >  0,  since  the  two  roots  are  symmetric  with  respect  to  the 
vertical  axis.  However,  the  RHS  of  Equation  5  does  have  two  roots  of  k,  the 
larger  one  turns  out  to  be  one  of  the  three  roots  for  the  negative  current 
case  (Figure  A-l).  If  mistakenly  chosen,  the  results  of  conversion  can  be 
very  different.  We  reasoned  that  since  a  contains  the  information  of  the 
direction  of  wave  propagation,  in  whatever  coordinate  system  chosen,  if  o  is 
for  the  forward  propagating  waves,  -a  is  for  the  backward  propagating 
counterpart.  For  a  reasonable  wind  condition,  the  major  propagation  direction 
of  short  waves  should  be  obvious.  Putting  equal  weights  on  +0  and  -0 
implies  equal  waves  travel  in  either  direction,  which  is  against  our  concept 
of  wind  wave  generation.  The  minus  sign  was  therefore  eliminated,  and  the 
directional  information  was  left  to  0. 

One  more  complication  for  the  reverse  current  case  (u  <  0)  should  be 
mentioned.  As  shown  in  Figure  A-l,  there  are  regions  (in  frequency)  with 
multiple  roots  of  wave  numbers.  Translated  into  the  matrix  equation,  this 
means  that  for  the  region  of  monotonic  k-a  relation,  the  submatrix  is 
square,  but  for  multiple  roots  region,  the  submatrix  is  not.  Augmentation  of 
this  rectangular  submatrix  is  necessary  to  enable  matrix  inversion.  This  can 
be  done  if  more  measurements  are  available.  For  example,  nx  and  measure¬ 
ments  provide  three  independent  set  of  spectra:  S^,  S^,  and  and 
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k  (radian /cm ) 


can  supply  enough  equations  for  as  many  as  three  roots  in  the  dispersion 

curves.  If  high  frequency  wave  height  measurement  is  also  available,  two  more 

independent  spectra,  can  be  added.  (S  .  S„  ,  and  S „  are  linearly 

r  r  ’  nnx  nny  n*  nx  ny 

dependent  when  wind  and  current  are  in  the  same  direction,  since  then 

A 

(S^  +  *  k^S^).  An  alternative  to  the  matrix  augmentation  is  to  select  a 

representative  wave  number  for  each  multi-root  frequency  (this  can  be  a 
weighted  wave  number  since  the  weighting  function  B  in  Equation  9  gives  the 
information  of  the  relative  contribution  of  each  root  to  the  total  spectral 
energy)  and  still  maintain  a  square  matrix  A.  The  alternative  method  is 
useful  when  not  enough  measurements  were  available  or  the  data  quality  of 
various  measurements  are  much  different,  and  propagation  of  error  into  the 
high  wave  number  region  is  not  desired.  The  latter  method  was  used  in  the 
field  data  conversion  presented  in  the  text.  The  following  shows  the 
conversion  of  the  hypothetical  spectrum  using  matrix  augmentation. 


Figures  A-2  to  A-5  show  the  pairs  of  spectra  for  u  ■  0,  40,  -60  and 

-100  cm/s.  In  each  figure,  the  left  hand  side  (plot  a)  are  the  frequency 

spectra  directly  calculated,  the  right  hand  side  (plot  b)  is  the  inverted  wave 

-4 

number  spectra.  The  recoverage  of  k  spectra  is  readily  seen. 


Figure  A3.  Converion  of  Wave  Frequency  Spectrum  (Left) 
to  Wave  Number  Spectra  (Right).  U  is  Current  Magnitude. 
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Figure  A4.  Conversion  of  wave  Frequency  Spectra  (Left) 
to  Wave  Number  Spectra  (Right).  U  is  Current  Magnitude. 


A5.  Conversion  of  Wave  Frequency  Spectra  (Left) 

:  Number  Spectra  (Right).  U  is  Current  Magnitude. 


APPENDIX  IV 


The  propagation  of  short  surface  waves 


on  longer  gravity  waves 


Michael  S.  Longuet-Higgins ,  Department  of  Applied 
Mathmatics  and  Theoretical  Physics,  Silver  Street. 

Cambridge,  England  and  Institute  of  Oceanographic 
Sciences,  Wormley,  Surrey. 

Abstract 

This  report  summarizes  some  theoretical  calculations 
carried  out  on  behalf  of  the  TOWARD  program  between  October 
1 9  8  A  and  September  1985.  The  purpose  was  to  determine  the 
theoretical  variations  in  the  wavelength  and  steepness  of 
short  gravity  waves  propagated  over  the  surface  of  a  train 
of  longer  gravity  waves  of  finite  amplitutde.  Such  varia¬ 
tions  may  be  calculated  once  the  orbital  accelerations  and 
surface  velocities  in  the  longer  waves  have  been  accurately 
determined  -  a  nontrivial  computational  task. 

The  results  show  that  the  linearised  theory  used  pre¬ 
viously  for  the  longer  waves  is  generally  inadequate.  The 
fully  nonlinear  theory  used  here  indicates  that  for  longer 
waves  having  a  steepness  parameter  AK  =  0.4,  for  example, 
the  short  wave  steepness  can  be  increased  at  the  crests  of 
the  longer  waves  by  a  factor  of  order  8,  compared  with  its 
value  at  the  mean  level.  (linear  theory  gives  a  factor  less 
than  2 )  . 

The  calculations  so  far  reported  are  for  free,  irro- 
t  a  t  Ion  a  1  gravity  waves  travelling  in  the  same  or  directly 
opposite  sense  to  the  Longer  waves.  However,  the  method 
of  calculation  could  he  extended  without  essential  diffi¬ 
culty  so  as  to  Include  effects  of  surface  tension,  energy 
dissipation  due  to  short-wave  b renting,  surface  wind-drift 
currents,  and  to  arbitrary  angles  of  wave  propagation. 


1.  Introduction 


An  important  component  of  radar  backscatter  from  the  sea  surface 

arises  from  the  Bragg  scattering.  This  involves  surface  wavelengths 

of  the  order  of  a  few  centimetres  for  X-band  radars,  or  tens  of 

centimeters  for  L-band .  In  both  cases  the  wavelengths  are  usually 

small  compared  to  the  dominant  wavelengths  of  ocean  surface  waves 
3 

(10  to  10  m) .  So  it  becomes  an  important  question  to  study  how  the 
short  wave  energy  is  distributed  with  respect  to  the  phase  of  the 
longer  waves 

In  the  present  study  we  shall  consider  the  classical  model  of  a 
short  train  of  gravity  waves,  of  small  but  variable  steepness  (\,k  , 
propagated  over  the  surface  of  a  longer  train  of  gravity  waves  of 
finite  steepness  AK.  as  in  Figure  1.  Early  workers  (Longuet-Iliggins 
and  Stewart,  1960)  assumed  that  A  K  I  ,  and  in  that  case  it  was 

found  that  the  variation  in  the  wavenumber  fa  and  amplitude  C.  of  the 
short  waves  was  given  by 

k/k 

a  jcx 

where  fa  and  a  are  the  (constant)  values  of  ft  and  CL  at  the  mean 
outface  level  and  ~y  =  j\  ^  ->  -  —  f  J  is  the  phase  of  I  In- 
long  waves.  This  gives 

^k/cik  ,  I  -i-^/jlv  /■  1  o(ak),  ,,.2) 

IV-  2 


-  /  -h  A  K  y-  -fa  C  (  a  k)  1 
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showing  that  the  short  waves  are  both  shorter  and  steeper  on  the  crests 

of  the  longer  waves  ! =  3.  h  .  However,  since  A  k  (  (AK)  =  0-4 -s 

'  V  /W»«5r 

the  maximum  steepening  predicted  is  less  than  2. 

Longuet-Higgins  and  Stewart  (1960)  interpreted  this  result  by 

assuming  (1)  that  the  phase  of  the  short  waves  was  conserved,  i.e.  that 


k<j,  - 


Or 


—  o 


V.<x 


t*.  v\  r 


(1.3) 


where  is  the  particle  speed  in  the  long  waves  as  seen  by  an 
observer  travelling  with  the  long-wave  speed  G ,  and  O'  is  the  intrinsic 
frequency  of  the  short  waves  in  a  frame  moving  with  speed  .  Next, 

(2)  that  the  intrinsic  frequency  O'  and  local  wavenumber  /?  of  the  short 
waves  were  related  by 


O-’.j'k 


a. a) 


where 


3 


was  the  effective  value  of  gravity  for  the  short  waves  i.e. 

(1.5) 


f  =  J  dV//*t 

W  being  the  vertical  component  of  orbital  velocity  in  the  long  waves; 
and  (3)  that  the  short-wave  energy  density  E  was  given  by 


£  *  j  a  +  4  **  aV  k  ,  (1>6) 

representing  the  potential  and  kinetic  energies  respectively.  The 
changes  in  shortwave  energy  El  over  the  long  wave  could  then  be 
attributed  to  (a)  advection  by  the  long-wave  orbital  velocities, 
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0\J 


together  witli  (b)  work  clone  by  the  straining  of  the  long  waves  against 
the  radiation  stress  of  the  short  waves. 

Garrett  (1967)  suggested  that  the  same  results  (1.1)  could  be 
interpreted  in  terms  of  the  conservation  of  wave  action 


where 


E 


/ 


(1.7) 


(1.8) 


is  an  alternative  form  of  the  short-wave  energy  density,  and  be 
introduced  the  equation 


9A 


3 


<  »- 


if?-',)'1] 


where  C  is  the  group-velocity  of  the  short  waves  (  C 


(1.9) 


i 

3 


c 


)■ 


Finally,  Bretherton  and  Garrett  (1968)  proved  the  validity  of 
equation  (1.9)  for  a  general  class  of  situations  where  a  group  of 
linearized  short  waves  of  wavenumber  k  is  propagated  through  a 
slowly  -  varying  medium  with  local  velocity  (j  ,  under 
the  general  assumption  that 

«  k),  If  «  a-f,  n.10) 

the  energy  density  E.  being  defined  as  if  the  medium  were  locally 
uniform. 

The  great  advantage  of  this  formulation  is  its  relative  simplicity, 
and  that  there  is  no  explicit  restriction  on  the  steepness  Alt  <* 
the  long  waves;  it  appears  necessary  to  assume  only  that 
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h  »  K  . 


d.ii) 


a .k  «  I  . 

) 

In  the  case  of  A  /Cf  inite ,  one  would  take  ns  the  effective  (vector)  gravity 


$ 


-  <V 


(1.12) 


where  e(  is  the  orbital  acceleration  in  the  long  wave. 

This  principle  has  been  partly  applied  (in  principle)  by  Phillips 

(1981)  to  calculate  the  variation  in  amplitude  of  short  capillary-gravity 

waves  riding  on  longer  gravity  waves.  The  calculation  could  not  be 

carried  through  in  detail  because  the  effective  gravity  was  not 

at  that  time  known  with  sufficiant  accuracy.  However,  the  accurate 

calculation  of  accelerations  in  steep  gravity  waves  has  recently  been 

carried  out  by  Longuet-Higgins  (1985),  and  from  this  it  is  possible 

to  infer  by  (1.12),  hence  both  the  shortening  and  steepening  of 

the  short  waves.  In  this  contribution  we  apply  the  results  to  short 

gravity  waves,  in  the  first  place,  with  application  particularly  to 

backscattering  in  L-band.  One  significant  result  is  that  for  finite 

values  of  A  K  the  short-wave  steepening  can  actually  be  much  greater 

than  that  given  by  linear  theory.  Moreover,  it  will  be  seen  that  the 

/ 

basic  calculation  of  Q  opens  the  way  to  the  solution  of  much  more 
realistic  cases,  including  the  effects  of  capillar!  ty,  short-wave  energy 
dissipation,  etc. 

The  assumptions  (1.11)  make  clear  that  the  principle  of  action 
conservation  is  only  approximate.  To  test  the  principle  in  a  simple 
case  we  study,  in  section  6  of  the  present  paper,  a  simple  model 
analogue  in  which  the  equations  of  motion  can  be  integrated  precisely 
by  numerical  methods.  The  degree  to  which  the  principle  is  satisfied 


is  reassuring. 
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2.  Formulation  of  the  problem 


Relatively  short  gravity  waves,  of  local  height  2<x  and  wavelength 
2i</k  .  ride  on  longer,  progressive  gravity  waves  of  finite  height 

and  wavelength  in  deep  water,  where  h»K  (see  Figure  1) . 

It  is  required  to  find  t\  and  O.k  as  functions  of  the  phase  of  the  long 
wave . 


The  intrinsic  frequency  cr  and  the  wavenumber  /?  of  the 

waves  are  assumed  to  be  related  by  equation  (1.4),  where  ^ 

/ 

magnitude  of  the  effective  acceleration  4  given  by  (1.12). 

is  always  normal  to  the  surface  of  the  longer  waves.  The 
O'  and  the  phase-speed 


short 
is  the 
Clearly 
f  requency 


(2.1) 


are  taken  as  positive  or  negative  according  as  the  short  waves  travel 
in  the  same  or  opposite  direction  to  the  long  waves.  denotes  the 
particle  speed  at  the  surface  of  the  longer  waves,  as  seen  in  a  frame 
of  reference  moving  with  the  long-wave  phase  speed  C.  In  this 
reference  frame  the  long  waves  appear  steady  and  the  free  surface  is 

,  we  have  <j  —  C  ;  (see  Lamb 

1932,  C.9) . 

To  determine  the  wavenumber  f\  at  points  along  the  surface  of  the 
long  waves  we  assume  that  the  phase  of  the  short  waves  is  conserved, 
that  is 


a  streamline.  At  the  mean  level  U  — 


*(%  -  0 


1  \  4  f  n  i ,  / 
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(2.2) 


say.  Hence 


c*  =  y/k  =  /;Yrc) 


or 


c  +  yj 


C  - 


rj  % 


'0=0 


(2.3) 


(2.4) 


a  quadratic  equation  for  C  with  solutions 

c  =  -/s  ±(ji\2py'\  A -»/;•  «-5> 

Having  found  C  we  may  calculate  /?  from  (2.3)  in  the  form  fa  ~  J  / C* . 

To  determine  the  local  wave  amplitude  (X  ,  we  assume  that  action 
is  conserved,  that  is  equation  (1.9).  In  the  steady  flow  relative 
to  the  moving  frame  of  reference  this  implies  that  the  flux  of  wave 
action  is  a  constant,  i.e. 

/ 


(<j,  -  Cj)  E  /  <T  “  U.J\:fa 


(2.6) 


where  C  the  group-velocity  of  the  short  waves  ( —  ^  c 
•  v 

and  E.  is  the  intrinsic  energy-density  of  the  short  waves,  given 
by  equation  (1.8).  Since  O'  =  |  equation  (2.8)  can  also  be 

written 


or 


(c.f,  Longuet-Higgins  and  Stewart  1961). 


(2.7) 


(2.8) 
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Finally,  the  local  wave  steepening  is  defined  as 


Jl  —  c\.  lx  !  (  cT  h  ) 


(2.9) 


where  a  bar  denotes  the 


values  at  the  mean  level 


Clearly  the  above  approach  depends  critically  on  the  accurate 
evaluation  of  the  velocity  and  the  orbital  acceleration  ftf  in  a 
(long)  gravity  wave  of  finite  amplitude. 


I 

I 
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3.  Method  o£  Calculation 


The  real ,  or  orbital  acceleration  in  a  surface  wave  must  be  care¬ 
fully  distinguished  from  the  apparent  accelerations  as  measured  by  a 
fixed  vertical  wave  gauge  (see  Longuet-Higgins  1985c).  The  real  accelera¬ 
tions,  both  vertical  and  horizontal,  vary  much  more  smoothly  than  the 
apparent  accelerations,  which  can  be  very  non-sinusoidal . 

Numerical  values  of  the  real  acceleration  were  calculated  by 
the  method  of  Longuet-Higgins  (1985a)  which  makes  use  of  a  set  of 
quadratic  relations  between  the  Four ier  coe f f icients  O.  in  Stokes's 

n 

series  for  the  cartesian  coordinates  in  terms  of  the  velocity 

potential.  Thus  if  ft  ~  ~  I*  t^ie  ^reG  surface  is  given  by 


<v> 


K 


J 


—  (X 


+  Z 


<x> 


<"  f  A ) 


(3.1) 


fa*  ~  //<•'  +  2  ah  (h^/c) 

where  <j>  is  the  velocity  potential,  then  the  coefficients  O~0  }  (  ^  CV  ^  ^ 

satisfy  the  relations 


cx 

o 

K 

+ 

cv.  s 

+  (rL 

+-  1  r 

‘r  +  -  ■  -  - 

< 

-  C 

0 

+ 

tl0  S 

1  r'i 

■*  ’l 

S'-  = 

0 

K 

+ 

+  S 

"f  ('<• 

1  6  J  -(■  -  -  -  = 

0 

with 

b 

h 

rr 

h  an , 

h  >  0  and  / 

v  1 

.  These  reations 

may  be 

(3.2) 


and  accurately  solved  for  a  given  value  of  the  phase-speed  C  (in  general) 


IV-10 


or  of  the  wove  amplitude 


A  - 


<x(  -h  i  \ 


(3.3) 


The  speed  q  at  tlie  free  surface  is  then  found  from  t lie  Bernoulli  relation 


(3.4) 


and  the  vector  acceleration  Y  from  the  general  relation 


2  mA' Mi  '  -}6*a  2Aa 

wher cjZ-X+ty  and  k  r  ^  H  y>,  the  stream  stream  function.  (A  star 


A  AA 


(3.5) 


denotes  the  complex  conjugate) .  In  real  terms  this  is 


*  "  ~f(Xt  +7/)VxW  -  (3.6, 


The  effective  gravity  O  is  then  found  from  (2.2). 

Because  of  the  slow  initial  rate  of  convergence  of  the  series  (3.1) 


at  high  values  of  a/y,  care  must  be  taken  to  include  enough  terms  in 


these  series.  A  recent  study  (Longuet-Higgins  1985b)  has  shown  tliat 


—  ^  j 

after  an  initial  rate  of  convergence  like  h  ,  CX  ultimately 

n 


converges  exponentially, the  transition  to  exponential  behaviour  occurs 


when 


<1  =•  »'c  =  O  (<r  JJ  »  where 


J-o  |/Ur  -  A/vj^J  .  <3-7> 


Since  individual  terms  in  the  differentiated  series  for  <X  ,  .  and  < 


at  first  increase  like  h  ,  it  is  important,  in  order  to  ensure 


sufficient  accuracy  in  the  calculation,  to  include  terms  with  a  some¬ 


what  in  excess  of  . 


Surface  profiles  corresponding  to  /\  f\  ="  O-  t  C)  •  2  f  0-$  O  •  y 
and  the  limiting  value  A  *"  O •  f\  ^2  are  sl'OWI1  *n  Flr.ure  2. 
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4.  Results;  variation  in  short-wave  length. 

Using  suffices  1  and  2  to  denote  values  at  the  long-wave  crest 
and  trough  respectively.  Figure  4  shows  the  relative  shortening  /A1 

at  the  crests  of  the  long  waves  as  compared  to  the  mean  surface  level, 
in  the  three  cases  when  fa  =  3.  }  10  and  IOO  .  Similarly  faa  f  fa 
shows  the  lengthening  in  the  long-wave  troughs.  For  values  of  A  K 
up  to  0.2  the  three  curves  corresponding  to  fa  *s  e?  ,  lO  and  f  0  0 
are  almost  indistinguishable,  and  even  when  A  the  re  is  little 

departure  from  the  representative  curve  fa  =  10  ,  when  k,/k  =3-0  and 
fa*  I  fa  0  •  ?  2.  •  Thus,  the  short-v;ave  length  varies  over  a  range  of 
about  3’s  to  1.  Tills  Is  for  C  >0  *  when  the  short  waves  travel 
in  the  same  sense  as  the  longer  waves.  Figure  5  shows  a  similar  plot 
when  C  <  0  ,  and  the  short  waves  travel  in  the  opposite  sense.  Here 
the  variation  in  fa  is  only  slightly  less.  However,  as  A /C 
it  can  be  shown  that  /?,  A  ->  t> o  when  C  >  0  *  but  remains 
finite  when  C  <  0  . 


I 


5 

8 
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waves,  as  compared  to  Lite  moan  level: 


As  Figure  5,  but  for 


l 


% 


V; 

v, 


5.  Variation  in  the  wove  steepness 


Figures  6  and  7  show  the  variation  in  steepening  of  the 
shorter  waves,  in  a  similar  manner  to  Figures  4  and  5.  Thus 

a,  -  -y  ,  •'>,  *  -  (5.D 

O.  R  CX  h 

Again  the  three  curves  corresponding  to  fe  •-=  SI  /Qand  /OO  lie 
very  close  together,  showing  that  not  only  the  wavelength  variation 
but  also  the  steepness  variation  is  practically  independent  of  short¬ 
wave  length. 

l^lien  however,  the  short-wave  steepness  may  vary 

by  a  factor  of  as  much  as  8  between  the  long-wave  crests  and  the 
mean  level.  This  compares  with  a  factor  less  than  2  given  by  linear 
theory . 

The  variation  of  steepness  Si.  over  the  profile  of  the  long  waves 
is  shown  in  Figure  8  as  a  function  of  y j\-.  ,  and  for  different 
values  of  AK.  using  the  representative  short  wavenumber  /?-$*. 
Comparing  AK  -  Q'1  with  AK-  0  •(  one  sees  the  distorting  effect 
of  nonlinearity  in  the  long  waves. 

Finally,  in  Figure  9  the  three  curves  of  Figure  8  are  plotted 
instead  against  »  tlm t  is  the  vertical  height  above  the 

mean  level  lj  .  It  now  appears  that  all  the  curves  collapse  almost 
onto  a  single  curve.  This  property  may  be  useful  in  approximate 
analytical  work.  The  appropriate  nonlinear  steepening  is  quite'  different 
from  the  linear  theory,  shown  in  Figure  9  by  the  broken  curve. 
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The  relative  steepening  oE  short  waves  at  the  crests 


1 


6.  A  pendulum  model 


As  a  check  on  the  principle  of  action  conservation,  which  was 
used  in  the  calculation  of  the  wave  steepness,  a  simple  mathmatleal 
model  was  constructed.  The  long  wave  was  represented  as  a  single 
"orbiting"  pendulum  of  length  /_  describing  complete  orbits  in  a 
vertical  plane  (Figure  10).  To  the  bob  of  the  longer  pendulum,  of 
mass  /'“'t ,  is  attached  a  smaller  pendulum  of  length  J!  and  mass  IY\  , 


where 


l  «  L 


n\  <<  f  | 


(6.1) 


It  is  convenient  to  choose  units  so  that 


L  =  / 


M  i  ».  -  I  , 


(6.2) 


then  the  kinetic  and  potential  energies  of  the  system  are  respectively 


T  -  J-  0  +  J  >vy  B  Gt.  >  ~  f )  -J-  J  y* 

~  (c~~>  &  +  ij  ~h  a \y\  f  c- o  ^  *f- 1  ^ 


(6.3) 


where  (?  and  <J)  denote  the  angles  which  the  longer  and  shorter  pendulums, 
respectively,  make  with  the  vertical.  The  Lagrange  equations: 
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f* 


^  T 

W 

(  •?*>  / 

2$ 

(  tT  \ 

IT 

2V 

l  ^  / 

> .() 

(6. A) 


yield 


§  +-  / in  [ •j'  '/  (Q  -  / )J  -  ■J"'  $ 

I 

Jj  +  [P  /*-*{)  ^  f , 


•  » 


These  can  be  solved  for  $  and^in  terms  of  t'^  ^  and  <f)  and 
integrated  numerically,  for  various  initial  conditions.  As  a  check 
on  the  time-stepping  we  may  monitor  the  total  energy 


T  +  V  -  H  . 


(6.6) 


When  either  i.  or  vanishes,  the  system  reduces  to  a  simple 
(long)  pendulum  of  length  L-/  obeying  the  equation  of  motion 


(6.7) 


J  $  -I-  (9  --  OlM j,  ^  -f-  /  ^  (6.8) 

where  is  the  value  of$at  the  highest  point  0 ,  <3  .  Writing  |v  «  /  /Q 


we  easily  obtain 
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(6.9) 


t  -f 


8 


/r  <(  0 


•/i  > 


°  fl  K 

a  known  elliptic  integral. 

/ 

It  will  be  seen  that  the  effective  gravity  is  always 

directed  radially,  and  has  magnitude 

j '  -  fj '  -  j  >  O  -  Q * ~  3  ^  >  (9 .  (6.io) 

_  a  t 

When  ^  is  always  directed  r a d  i  a  1 1  y  outwards . 

We  are  interested  in  conditions  when  O I  so  that  the  long 
pendulum  describes  complete  orbits  rapidly,  and  gravity  produces  a 
relatively  small  perturbation.  Then  I  »  and  expansion  of  (6.9) 

in  powers  of  It  yields 

Q(~  =r  -h  0  +  Offs*)  (6.H) 


and  hence 


S  -  erfc  —  (x.  q  [  f~)  {lx  ) 


(6.12) 


where 


Cr  ~ 


o 


8.. 


/-/r1 


/  -  K 


(6.13) 


a  I  - *  . 

Equation  (6.12)  represents  a  uniform  rotation  together  with  a  periodic 
perturbation,  of  radian  frequency  c  given  by  (6.13).  In  the  next 
approximation  we  find 


9 

O 

i  *  7  *i 

I  -  Hr  I  ,  ir 

'/• 


(6.16) 
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In  the  general  cnse  when  £  >  O  and /M  >  CJ  we  may  define  the 


energy  of  the  small  pendulum  as 


E.  -  j  ^  ^  -f-  £  in  I r'  y  ( (/«  i  /\)  +  I J 


(6.15) 


where 


is  the  effective  gravity.  Then  using  equation  (6.5)  and 


neglecting  terms  in  7*1  we  find 


=  J.m  1^3  &  4  m  $  -  J  (9  -cJm-.  r’  (jb  -  0) 


-  8(9-  «8)  //  -Vj] 


(6.16) 


It  will  he  noticed  that  tiie  terms  in  ^  arc  absent.  After  averaging 
over  the  fast  phase  ^  we  obtain,  to  lowest  order 


/  c*  j 

r:/r 


>  =  J  -<V  <9 


•  n 

V'»-  (, 


(6.17) 


as  the  approximate  variation  in  the  short-wave  enprgy. 


An  alternative  definition  of  the  short-wave  energy  is  to  replace 


the  right  hand  side  of  (6.15)  by  (‘f'-O')  »  so  taking 


account  of  the  rotation  of  the  longer  pendulum.  We  denote  the  corre- 


EU 

.  This  leads  to  an  equation  similar  to  (6.17) 


To  test  these  two  definitions  of  the  energy  we  investigated  the>fcase 


1=10 


(6.18) 


with  initial  conditions 


8  -  f  -  o 


9  -  l 


ft  -  VC 


(6.19) 


at  £*  =  0  •  Figure  11  shows  the  trajectory  of  the  mass  In,  and  Figure  12 

»  •  • 

indicates  the  angular  velocity  '/  as  a  function  of  $  ;  the  point  &  ~  ( 
corresponds  to  -r  Q  ,  and  B  =  £-2 /S’  tP  d  “  1\  (lowest  point  of 
the  larger  mass  /VJ).  Clearly  the  f  Iiic  I  n;U  ionsin  (j)  are  as  great 
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as  18%.  Figure  13  (upper  curves)  show  the  limits  in  the  variation  of 


Ei  ,  which  are  about  7  percent  of  the  mean.  The  lower  curves,  on  the 


r-  "  ■ 

other  hand,  show  the  fluctuations  in  LZ  ,  which  are  clearly  much  less  — 

rr  /f 

about  0.5  percent.  Thus  £_  is  a  more  satisfactory  definition  of  the 


energy  of  the  smaller  pendulum  than  is  IE.  .  However  the  mean  value 


rt  ■ 

)f  C.  still  varies  by  about  20%  over  a  complete  cycle  of  L7. 


Finally,  to  test  the  action  priniciple  we  calculated  the  two 


cr  '/  /  tr'7  0  '  " 

itios  CZ  / H  and  C  f  (5  where  (J  and  <7  are  given  by  expressions  analogous 


to  (6.14),  viz 


cr  r  - 


/I  ’ 


(6.20) 


a"-  - ^ - 4  ,  <?'-  (4-9)  *  •  I/O  ■ 

I  -  ly,,7+ 1  ft"*  )  9-f  =  r  i  (6.21 

r-  "  /  " 

The  results  are  shown  in  Figure  14.  The  total  variation  in  tL  /  rj  is 
clearly  less  than  that  of  E /(J  and  is  of  the  order  of  4%  only.  Some 


tx  =  I/Q 


(6.21) 


of  this  variation  may  be  due  to  the  approximate  nature  of  the  frequency 


formula  (6.14) 


This  example  was  typical  of  those  that  were  tested.  On  the  whole. 


therefore,  the  print:  i  pi  eof  action  conservation  was  well  verified. 
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Figure  1 1 


Trajectory  of  the  mass  vhen  -  ft)"' 

_  /f  ’ 

/h  r  (Q  ^  and  the  Initial  conditions  are  that 

(9  -  <f>  -  O.  0^1,  ji  *•  O  when  C~  -  O  . 


IV-29 


Maximum  and  minimum  values  of  the  energy 
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Figure  I.  Definition  sketch  for  short  waves  on  long  waves, 


Figure  2.  Surface  profiles  of  gravity  waves  in  deep  water. 
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Figure  4.  The  relative  shortening  of  short  waves  at  the  crests 
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Figure  5.  As  Figure  3,  but  for  t.  <  o  . 


Figure  f>.  The  relative  steepening  of  short  waves  at  the  crests 


(  )  and  in  the  troughs  (-1  l  )  of  long  waves, 


as  compared  to  the  mean  level:  C  >  O  . 


Figure  7.  As  Figure  6,  but  for  C  < 


Figure  8.  The  relative  steepening  ^  as  a  function  of  tbe 


horizontal 
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Figure  9.  The  relative  steepening  as  a  function  of  t  lie 


vertical  coordinate  —  y  J  j /.  when  I V  y ,  c  <  O 


Figure  ID.  Definition  sketch  for  the  pendulum  model. 
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1 .  In  t  r oduc  t ion 


Meteorological  measurements  at  the  NOSC  tower  were  collected  in 
order  to  obtain  accurate  estimates  of  the  wind  stress  at  the  air- 
sea  interface.  Two  independent  systems  were  used  to  measure  the 
stress:  sonic  anemometer,  which  directly  measures  the  momentum 
flux,  and  hot  film  system,  which  indirectly  evaluates  the 
momentum  flux  based  on  a  turbulent  energy  balance.  Because  the 
flux  data  presented  in  this  report  are  based  on  the  hot  film 
system,  the  theory  required  to  obtain  the  wind  stress  from  this 
method  will  briefly  be  described.  Data  from  the  sonic  anemometer 
are  not  yet  ready  for  distribution. 

The  basic  meteorology  quantities  measured  were  wind  speed  and 
direction,  air  temperature  and  dew  point  temperature.  The  water 
temperature,  meausured  by  NOSC,  was  incorporated  into  the  meteoro¬ 
logy  data  set.  In  addition  to  the  momentum  flux,  high  frequency 
measurements  at  temperature  and  water  vapor  were  collected  in 
order  to  evaluate  the  fluxes  of  latent  and  sensible  heat.  During 
TOWARD  85,  the  meteorology  data  set  was  supplemented  by  twice 
daily  radiosonde  profile  data,  launched  from  the  tower.  Net 
solar  and  infra-red  radiation  measurements  were  additionally 
collected  during  spring  1985.  The  measurements  are  summarized 
in  Table  1. 

Estimates  of  the  momentum  flux  by  the  dissipation  method  (hot 
film)  were  continuous  during  normal  tower  operations.  Sonic 
anemometer  measurements  were  restricted  to  conditions  such  as 
moderate  to  high  westerly  winds  and/or  times  of  overflights. 

The  evaluation  of  data  quality  required  consideration  of  mast  flow 
distortion,  fog/rain  conditions,  and  convection  off  the  tower  body 
during  light  winds.  Because  the  mast  on  which  all  the  meteorlogy 
sensors  were  mounted  (except  radiation)  was  located  on  the  north 
end  of  the  NOSC  tower,  the  greatest  flow  distortion  occurred 
during  southerly  flow.  The  most  favorable  direction  for  measure¬ 
ments  were  northerly  winds  from  the  west  through  the  northeast. 
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Fog  and  rain  interfered  considerably  with  both  the  sonic  and  hot 
film  systems.  Because  the  sonic  is  based  on  the  relative  time 
for  a  sound  pulse  to  travel  over  a  short  distance,  a  thin  film 
of  water  on  either  the  source  or  detector  gave  an  erroneous 
measure  of  the  wind  velocity.  The  hot  film,  on  the  other  hand, 
when  in  contact  with  a  water  droplet,  produced  a  voltage  spike 
with  the  latent  heat  release  during  evaporation. 

Finally,  light  winds  adversely  affected  the  quality  of  the  flux 
estimates  due  to  both  enhanced  flow  distortion  and  convection  off 
the  body  of  the  tower. 

Tabe  1  1 


Arrangement  of  sensors  on  the  NOSC  Tower 
(variables  defined  in  Section  2) 


Sensor  Variable  Measured 


Sonic  anemometer 
Lyman-a 
Cold  wire 
Hot  film 

Tair>  Tdp 

Cup  anemometer 

Wind  vane 

Water  temperature 

Solar / Inf rared  Rad 

Radiosondes 


u.v.w.T 

q  »q ' 

T« 

e 

Ta ir  »  ^ 

u  (slow  response) 
9  (0°  -  North) 

To 


Height  above  MSL 
22  m 
19  m 
22  m 
2  2  m 
18  m 
18  m 
17  m 
-3  m 


2 .  Basic  Theory 

2.1.  The  direct  method 

Generally,  the  wind  stress  (or  momentum  flux),  T,  may  be  written 
in  nearly  exact  terms  as: 

t  *  -  puw  (  1  ) 


where  p  is  air  density,  and  u  and  w  are  the  instantaneous  horizon¬ 
tal  and  vertical  wind  velocities.  If  an  averaging  time  is  chosen 
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to  be  long  enough  so  that  w  =  0,  then  equation  (1)  reduces  to 


i  =  -  p>u  1  w' 


(2) 


The  sonic  anemometer  uses  this  method  for  obtaining  The  fric¬ 
tion  velocity  is  simply  a  scaling  parameter  for  wind  stress  (with 
units  of  velocity)  and  is  defined  as: 

I 

u*  =  (  x/p)  ’/2  (3) 


Because  geophysical  modellers  of  air-sea  interaction  processes 
needed  spatial  field  estimates  of  the  wind  stress,  the  bulk  aero¬ 
dynamic  formula  was  derived.  This  formula  requires  knowledge  of 
the  mean  meteorological  quantities  (winds  and  temperatures)  and 
is  written  as : 


t 


PCDu2 


(A) 


where  CD  is  the  drag  coefficient.  The  drag  coefficient  is  known 
to  vary  considerably  with  atmospheric  stratification  and  height. 

If  the  measurement  height  is  chosen  to  be  10  m  above  the  surface 
and  the  stratification  is  neutral,  C3  'v  .001,  i.e.,  roughly  a 
constant  up  to  a  windspeed  of  10  m/sec.  Above  10  m/sec,  Cp 
becomes  a  slowly  increasing  function  of  windspeed.  As  the  atmos¬ 
pheric  surface  layer  becomes  more  unstable,  Cp  increases;  and 
conversely,  as  the  surface  layer  becomes  more  stable,  decreases. 

Wind  stress  data  are  often  reported  in  terms  of  the  drag  coeffi¬ 
cient  and  its  neu t r a  1- s t r a t i f ic a t ion  counterpart,  CDN.  For  given 
meteorol  igical  conditions,  Cpi;  may  be  evaluated  from  CD  and 
knowledge  of  the  surface  layer  stability,  assuming  that  the 
observed  windspeed  is  common  to  the  surface  layer  equations.  The 
only  reason  CDN  is  calculated  from  TOWARD  data  is  to  compare  to 
past  investigations  of  the  drag  coefficient  over  the  sea.  Even 
further,  since  stability  is  a  known  parameter  that  affects  the 
magnitude  of  the  drag  coefficient,  performing  statistics  on 
the  neutral  drag  coefficient  hopefully  will  yield  information 
on  a  secondary  parameter  that  affects  Cp,  i.e.,  sea  state.  Such 
an  analysis  is  more  simplified  by  stating: 


CD( *rWS)  =  CDN(WS) 


(5) 


where  t  ts  a  stability  function,  which  in  known,  and  WS  is  wave 

state,  which  is  unknown.  Any  trends  in  C_.„  with  wave  conditions 

UN 

will  provide  clues  to  this  function;  preliminary  results  are 
presented  herein.  Both  the  sonic  and  hot  film  systems  are  capable 
of  providing  estimates  of  U*  ,  Cp,  and  . 


2.2.  The  indirect  dissipation  method 


The  hot  film  system  utilized  the  dissipation  technique  for 
estimating  the  wind  stress.  This  method  requires  estimates  for 
the  balance  of  turbulent  kinetic  energy  (TKE) ,  where  the  TKE 
equation  is  written  as: 
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oz  Y 
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(6) 


where  e  is  TKE,  g  is  gravity,  e  is  the  TKE  dissipation  rate,  and 
p  is  pressure.  If  the  flow  is  not  excessively  stable  or  unstable, 
the  third  term  on  the  r.h.s.  in  (6)  is  relatively  small  and  can 
be  ignored.  Similarly,  for  steady  state  conditions. 


and  (6)  reduces  to: 
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If  one  assumes  that  the  vertical  profile  of  horizontal  windspeed, 
au/at .  may  be  defined  as: 

3u  u* 

-  =  -  tnU/L)  ,  (8) 
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where  ♦m  is  a  known  stability  function  dependent  on  z/L,  i.e. 
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then  the  friction  velocity  is  evaluated, 
to  be  : 

ekz 


♦h(EH 


after  rearranging  (7), 


(10) 


L  is  the  Monin-Obukhov  length,  k  is  the  von  Karman  constant. 

The  dissipation  rate,  e,  is  evaluated  by  spectrally  decomposing 
the  horizontal  wind  field  according  to  turbulent  eddy  frequency. 
The  inertial  subrange,  i.e.,  the  turbulence  with  frequencies 
of  motion  above  around  1  Hz,  is  well  described  by  the  Kolmogoroff 
spectrum ;  i.e., 

S(kx)  =  at2/3k^5/3 

(ID 

where  kA  is  wavenumber,  S(kA)  is  energy  density,  and  a  is 
a  constant  equal  to  0.52.  Equation  (11)  is  transformed  simply 
into  frequency  space  by  applying  Taylor's  hypothesis,  which 
assumes  the  concept  of  frozen  turbulence  advection,  i.e. 


2  nt 

k\  = - 

u 


Equations  (11)  and  (12)  combine  to  become: 


(  12) 


S  (  f  )  =  a(^-V5/3e2/3f-5/3 
u 


(13) 
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Although  the  spectral  decomposition  of  the  horizontal  windspeed 
will  yield  s(f)  and  f,  and  consequently  e,  the  method  used  in 
TOWARD  was  to  integrate  (11)  and  then  apply  (12).  This  yielded: 


a 


2 
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3or  cu  2/3 

2  2  it 


(C-2/3-c-2/3) 


(14) 


where  the  data  aquisition  system  employed  filter  frequency 
bounds,  fj  =  5  Hz  and  f^  =  50  Hz.  o^  is  a  quantity  evaluated 
due  to  the  filtering  process. 

With  (14),  e  was  evaluated,  and  consequently  with  (10),  estimates 

of  u*  were  provided.  These  estimates  are  provided  on  both  5  and 

30  minute  intervals  during  the  TOWARD  experiment.  For  those 

occasions  when  the  data  were  of  bad  quality,  u*  estimates  were 

_  3 

made  by  assuming  a  neutral  drag  coefficient,  Cpfj,  of  1.4x10 
(according  to  Large  &  Pond,  1981).  The  drag  coefficient,  CQ , 
applied  to  (4)  in  this  case  used  the  relation 

CD  =  (CDN~,/2  +  Vk)'2  H5) 


where  H1  is  a  stability  function. 
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>  C„„  and  for  stable  conditions 
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ties-  CD  *  CDN 


For  unstable  conditions 

C ,,  *"  C  ,  .  For  neutral 
U  UN 
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3.  Preliminary  Scientific  Results 


Plotted  in  Figure  1  are  the  dissipation-derived  neutral  drag 
coefficients  from  both  phases  of  the  TOWARD  experiment.  The 
conditions  of  these  data  include  westerly  and  northwesterly 
directions;  these  directions  were  dominant  in  the  data  set  and 
were  estimated  to  induce  minimal  flow  distortion. 

Figure  1  was  statistically  analyzed  at  1  m/sec  intervals;  these 
results  are  listed  in  Table  2.  Tabulated  alongside  the  results 
of  Figure  1  are  parallel  results  from  measurements  for  short 

fetch  (less  than  4  km). 
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The  results  of  Table  1  may  be  summarized  as  follows:  1)  Cn 

DN 

constant  over  the  windspeed  range  of  4  to  9  m/sec;  2)  The  scatter 

in  the  data  decreases  with  increasing  windspeeds;  3)  C:  .  increases 

DN 

substantially  as  a  function  of  windspeeds  for  winds  above  9  m/sec; 
and  4)  C  is  generally  larger  for  conditions  of  short  fetch  than 
for  long  fetch. 

The  data  analysis  indicated  systematic  trends  within  the  range 
of  windspeeds  where  was  found  to  be  constant,  i.e.,  4<U1Q<9 

m/sec.  CQN  was  observed  to  decrease  with  time  during  the  course 
of  a  sea  breeze  (see  illustration  -  Figure  2),  and  to  increase 
markedly  prior  to  a  frontal  passage  (see  illustration  -  Fiugure  3). 
Both  of  these  cases  are  assumed  to  be  due  to  the  surface  wave 

Table  2 


CDN  vs*  u10:  statistics 


Windspeed 

Long 

Fetch 

Short 

Fetch 

(m/sec)  Interval 

'°3  CUN 

°lCDtl ) 

io3  cDN 

o(ClJfJ) 

1-2 

2.12 

— 

5.03 

— 

2-3 

1.49 

0.83 

2.02 

0.85 

3-4 

1.43 

0.85 

1.66 

0.4  1 

4-5 

0.96 

0.44 

5-6 

0.98 

0.44 

6-7 

0.90 

0.38 

7-8 

0.86 

0.46 

8-9 

0.94 

0.25 

9-10 

1.41 

0.  10 

10-11 

1.48 

- 

From  (4-9)  m/sec;  103  = 

0 . 9  3  ±. 4 0 

field  undergoing  rapidly  changing  conditions.  In  the  sea  breeze 

case,  the  wave  field  is  rapidly  growing  at  all  wave  frequencies 
and  is  making  an  asymptotic  adjustment  towards  equilibrium  con¬ 
ditions.  In  the  second  case,  the  convergence  of  wave  trains  in 

the  pre-frontal  region  (from  both  ahead  and  behind  a  cold  front) 
acts  to  increase  the  steepness  of  waves  of  all  frequencies, 
thereby  increasing  the  surface  roughness.  Both  these  findings 
from  TOWARD  will  significantly  improve  our  understanding  of 
sea t t e r ome te r  backscatter  algorithms  and  their  relations  to  the 
overlying  winds. 

4.  Illustration:  31  October  1984 

In  Figure  4  and  5,  one  may  find  plots  of  meteorology  data  during 
the  SAR  flight  on  31  October  1985.  Because  the  flight  was  con¬ 
ducted  after  1300,  the  sea  breeze  was  westerly  with  windspeeds 
strong  enough  to  exhibit  confidence  in  the  data  quality. 


5.  Meteorology  during  SAR  flights 


Tabulated  in  Table  3  are  meteorological  data  quality  assessments 
and  general  conditions  during  each  SAR  flight. 

Table  3 


Flight 

2  Oct 

Wind 

Meteorology  during  over  Flights 
(*  =  CV  990) 

Conditions  Str a£ i E ication 

Di ss ipat ion 
Data  Quality 

S 

3-4 

N 

0 

6  Oct 

NW 

8-5 

S 

G 

7  Oct 

NW 

5-7 

N 

G 

17  Oct* 

NW 

7-8 

II 

G 

31  Oct* 

W 

3-6 

N 

G 

4  Nov* 

NW 

2-4 

U 

G 

7  Nov* 

W 

4-5 

S 

G 

13  Mar 

W 

4-5 

N 

G 

18  Mar 

NW 

3-10 

tl 

G 

19  Mar* 

W 

5.-6 

N 

G 

2  0  Mar 

NW 

5-7 

II 

G 

24  Mar 

W 

3-4 

S 

G 

2  5  Mar 

SW 

6-7 

S 

G 

G 

2  6  Mar 

SW 

4-5 

s 

27  Mar* 

SW 

5-8 

s 

I 
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South  West 

North  West 

Sou  th 

West. 

Windspeed  (m/sec) 

Stable 
lloiit  r  a  1 
Unstable 

V- 9 

Good 

Intermi ttent 

6.  Research  Plans 


6.1.  Short  Term 

a.  Technical  report  of  bulk  data  and  u*  estimates  from  the  dissi¬ 
pation  method.  Jan.  1986 

b.  Presentation  of  these  TOWARD  findings  described  herein  at  the 
Air-sea  Interaction  Conference,  Miami,  Florida.  Jan.  1986. 

c.  Submit  to  JGR  results  of  (b).  Spring  1986. 

d.  Cooperation  with  NPS,  NRL ,  UK,  JPL,  Riso  and  NOSC.  Ongoing. 

6.2.  Long  Term 

a.  Sonic  data  analysis  of  u* .  Spring  1986. 

2 

b.  Spectra  of  U,  0,  T,  g,  u* ,  w'T',  w'g'.  Spring  1986. 
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1 .  Background 

In  this  section  a  brief  overview  of  the  SAR  experiments 
and  theories  for  imaging  ocean  waves  is  presented.  The  exper¬ 
imental  results  and  the  limitations  imposed  on  them  are  reviewed. 
Predictions  of  the  various  SAR  imaging  theories  are  discussed. 

An  overview  of  SAR  flight  data  for  the  TOWARD  experiment  is  also 
presented. 

1.1.  Introduction: 

Images  of  ocean  waves  produced  from  the  MARINELAND  and 
MARSEN  experiments  have  shown  improved  contrast  or  Peak-to  - 
Baseline  Ratios,  PBR ,  by  adjusting  the  focus  parameters  in  the 
SAR  optical  processor.  Only  optical  recording  and  optical 
processing  were  used  in  these  experiments.  These  have  limited 
dynamic  range,  and  are  less  precise  in  determining  ocean  wave 
parameters  and  focus  settings  compared  to  digital  processing. 

It  was  considered  essential  in  TOWARD  84/85  that  digital 
processing  provide  the  basis  for  careful  assestment  of  the 
conflicting  SAR  imaging  theories.  Hence,  SAR  measurements  were 
recorded  in  both  optical  and  digital  modes.  Also,  the  images 
were  processed  in  both  optical  and  digital  modes. 

1.2.  SAR  FLIGHTS  in  the  TOWARD  EXPERIMENT 

The  CV-990  SAR  flight  pattern  over  the  NOSC  tower  are 
shown  in  Figure  1  which  shows  the  approximate  bearings  for  the 
legs.  The  flight  dates  and  duration  over  the  tower  in  phase-1 
are  given  in  Table  1.  The  environmental  conditions  during  these 
flights  are  shown  in  Table  2.  A  typical  CV-990  flight  pal  tern 
in  phase-1  is  outlined  in  Figure  2  and  the  JPL-SAR  characteristics 


CV-990  SAR  Flight  Paths  Over  Tower 


Table  1 .  Aircraft  Flights  During  Phase-I  of  TOWARD  84/85 
The  Flight  Intervals  are  Noted  for  Each 


CV-990 

F-4 

2  October  1984 

1046  -  1156 

6  October  1984 

1505  -  1600 

7  October  1984 

904  -  958 

17  October  1984 

1253  -  1536 

31  October  1984 

1257  -  1542 

1500  -  1614 

4  November  1984 

1512  -  1700 

1131  -  1231 

7  November  1984 

1356  -  1640 

Table  2 .  Environmental  Conditions  During  Flights  in  October- 
November,  1984:  Phase-I. 


Date,  1984 

17  Oct 

31  Oct 

4  Nov 

7  Nov 

Time  Analyzed 

1250- 

1300- 

1512- 

1412- 

1413 

1600 

1625 

1625 

Wind  Speed  (knots) 

15 

6-8 

5-7 

10 

Wind  Direction  (deg) 

~275 

~280 

~300 

~250 

Mean  Current  (cm/sec) 

<5 

<5 

<5 

~10 

(5  cm  Below  Surface) 

Significant  Wave 

126 

124 

158 

110 

Height  (cm) 

Air  T  (#C) 

17.2 

17.5 

18.5 

17.5 

are  given  in  Table  3. 


Ten  SAR  frames  have  been  selected  by  the  SAR  Data  Analysis 
Committee  for  focusing  studies.  These  are  shown  in  Table  4. 

The  frame  that  has  been  analized  is  for  leg  1  at  38,000  feet 
altitude,  collected  on  October  31,  1984.  The  processing  for  the 
other  frames  is  now  completed,  and  detailed  analysis  is  presently 
in  progress. 


1 . 3  Theoretical  and  Experimental  Background  on  Focusing  Effects 


The  available  theoretical  and  experimental  investigations 
on  focusing  is  reviewed  here  to  provide  the  perspective  for  the 
focusing  analysis  discussed  in  this  study.  Defocusing  has  been 
related  to  motion  of  the  ocean  surface.  Two  such  motions  are 
dominant:  1)  the  phase  velocity  of  dominant  ocean  surface  waves, 

and  2)  their  orbital  velocity.  According  to  Jain  (1978,  1981): 
the  optimum  image  is  obtained  when  the  matched-filter  of  radar 
processor  is  adjusted  in  proportion  to  the  azimuth  component  of 
the  wave  phase  velocity.  Also,  Harger  (1985  a,b)  refers  to  the 
surface-wave  phase  velocity  as  the  source  for  defocusing  of  the 
SAR  images,  and  hence  for  the  modeling  and  analysis  of  surface 
waves . 


Alpers  and  Rufenach  (1979)  state  that:  the  wave  like 
patterns  in  the  SAR  image  is  caused  by  orbital  motions  and  will 
occur  even  for  uniform  radar  scattering  cross  sections .  The 
motion  causes  degradation  in  azimuthal  resolution  and  is  dominant 
for  azimuthally  traveling  waves.  Also,  they  state  that  velocity 
bunching  can  cause  an  enhancement  of  ocean  wave  patterns,  or  the 
occurrence  of  two  and  possibly  three  intensity  peaks  per  true 
wave  length,  depending  on  ocean  and  SAR  parameters.  Their  work 
indicates  that  images  can  not  be  improved  by  refocusing  the 
optical  processor  of  the  radar.  Plant  and  Keller  (1983)  use  a 
two-scale  Radar  Wave  Probe  simulation  of  the  SAR  to  infer  that 


VI-5 


Table  3.  Nominal  SAR  Data  Processing  Characteristics  for  JPL  L-Band  SAR 


4  OPTICAL  Channels  -  HH  -  VV  -  HV  -  VH 
4  DIGITAL  SUBFRAMES  -  HH  -  W  -  HV  -  VH 
Angle  of  Incidence  at  Center  of  Image  ■  35° 

Optical/Digital  Slant  Range  Coverage  ■  7  km  (Quad  Pol  Mode) 

Optical/Digital  Slant  Range  Coverage  “  12  km  (Dual  Pol  Mode) 

Optical  Along  Track  Coverage  ■  60  km 

Digital  Along  Track  Coverage  (Raw)  ■  30  km 

Digital  Along  Track  Coverage  (Image)  **  11  km 

Digital  Image  Size  m  927  pixels  range  x  1024  pixelB  azimuth 

Optical  Image  Size  **  50  mm  cross  track,  1:250,000  along  track 

RFP  =  3  x  (CV-990  Ground  Speed  in  m/s) 

=  600  @  200  m/s,  750  0  250  m/s 

Digital  Image  Pixel  Spacings  *  7.5  m  Slant  Range  x  11.0  azimuth 


Table  4.  SAR  Frames  Especially  Processed  for  Focusing  Studies 


Days - 

- -17  Oct 

31  Oct 

4  Nov 

7  Nov 

Alt  (feet) 

7,000 

20,000 

Leg  2 

Legs  4,  7 

Leg  1 

Leg  2 

Leg  2 

38,000/38,000 

Legs  1,  2 

4,  7 

SAR  images  are  not  faithful  reproductions  of  cross-section 
variations  over  the  surface.  Distortions  in  the  image  are  due 
to  the  velocity  bunching  mechanism.  Application  of  focusing 
correction  does  not  remove  the  distortions  from  the  simulated 
image . 

Based  on  experimental  results,  Shuchman  and  Shemdin  (1983) 
indicate  that,  based  on  analyis  of  modulation  depths  as  a  function 
of  focus  setting  in  the  SAR  processor,  the  maximum  in  modulation 
depth  occurs  when  the  motion  correction  corresponds  to  a  velocity 
that  is  of  the  order  of  the  wave  phase  velocity  of  the  dominant 
surface  waves.  This  occurs  for  both  azimuth-  and  r ange - t r ave 1 ing 
waves.  Therefore,  the  detectability  of  these  waves  can  be 
improved  by  focusing.  Jain  and  Shemdin  (1983)  investigate  this 
point  in  the  MARSEN  data  using  both  radar  and  in-situ  measurements. 
They  show  that  the  focus  dependance  for  optimum  imaging  is  that 
for  a  surface  moving  with  the  wave  phase  velocity.  Also, 
azimuthally  traveling  waves  are  as  visible  as  range  traveling 
waves  when  the  optimum  focus  setting  is  made. 


2 .  Data  Processing 


In  this  section  analysis  of  the  digitally  processed  image, 
for  leg  1,  38,000  feet  altitude  on  31  October  1984,  provided 
by  JPL  is  discussed.  Examination  of  the  visual  images  at 
different  focus  settings  is  described.  The  procedure  for  locating 
the  data  blocks  on  the  visual  images  and  then  on  the  computer 
compatible  (cc)  magnetic  tapes  is  presented.  Definitions  of  the 
terms  and  parameters  used,  and  the  significance  of  the  information 
collected  are  discussed.  Processing  methods  that  are  not  optimal 
(in  the  SNR  sense)  are  compared  to  that  which  is  optimal. 

2.1.  Examination  Photographic  Representations  of  Digitally 
Procejssejj  Images 

A  sample  of  digitally-processed  image  provided  by  JPL  is 
shown  in  Figure  3  for  one  focus  setting.  JPL  provided  13 
such  images  at  focus  setting  ranging  from  -1.75c  to  +1.75c,  where 
c  is  the  wave  phase  velocity  of  dominant  surface  wave.  The  image 
in  Figure  3  is  therefore  referred  to  as  the  Black  and  White  (B&W) 
image,  for  later  comparison  with  computer-enhanced  results. 

i 

2.2  Locating  Data  Blocks 

Two  blocks  of  data  are  shown  in  Figure  ^  ,  each  block  is 
256  pixels  (10.98m  per  pixel)  by  256  pixels  (10.98m  per  pixel). 

One  block  (No.  2)  is  in  an  area  where  the  water  depth  is 
approximately  30  fathoms  (60m).  The  second  block  (No.  3)  is  in 
an  area  where  the  water  depth  is  approximately  13  fathoms  (26m). 
The  second  block  includes  the  NOSC  tower  where  the  wave  measure¬ 
ments  were  carried  out.  All  blocks  have  the  same  inner  and  outer 
ranges.  The  location  of  the  data  blocks  are  chosen  such  that 
they  will  contain  the  most  visible  waves,  as  shown  is  Figure  3. 

The  NOSC  tower  location  is  shown  in  Figure  5. 
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3.  Sample  of  Digitally  Processed  Image  Used 
in  Present  Analysis. 


2.3.  Locating  Data  Blocks  In  Data  Files 

The  tape  format  is  such  that  each  focus  file  has  1,280 
blocks  and  each  block  contains  1,024  bytes.  The  1,024  bytes 
compose  one  line  along  the  azimuth  (i.e.  constant  range)  direction. 
There  are  1,280  lines  parallel  to  the  azimuth  direction,  each  is 
at  a  unique  range.  The  lines  are  divided  into  4  equal  segments 
each  in  one  block.  Therefore,  each  line  has  256  bytes  in  a 
block.  Also,  each  block  has  256  lines.  The  distance  from  the 
edge  of  the  image  is  measured  to  the  edge  of  the  data  block.  This 
is  converted  to  number  of  pixels  using  the  scale  1  pixel  *  10.98m. 

The  number  of  pixels  obtained  is  the  number  of  lines  (or  blocks) 
to  be  skipped  before  reading  the  tape.  For  recording  the  data 
for  block  2,  the  first  256  points  of  each  line  within  the  range  limits 
of  the  block  are  skipped  and  the  second  256  points  (bytes)  are 
recorded.  In  recording  the  third  block  data,  the  first  512  points 
of  each  range  line  is  skipped  and  the  following  256  points  are 
recorded.  This  is  done  for  256  lines.  Same  procedure  is  used  on 
all  the  focus  files. 

2.4.  Averaging  Methods 

To  investigate  the  spectra  of  dominant  wave  numbers,  different 
smoothing  and  averaging  methods  are  used.  Space  averaging  refers 
to  averaging  values  of  a  few  pixels  in  the  range  direction  at  a 
given  range  value.  This  is  done  such  that  an  equal  number  of 
pixels  are  taken  on  each  side  of  the  given  range  number  ,  which 
is  taken  to  be  the  range  value  of  the  cut.  All  cuts  in  the  processing 
are  taken  to  be  parallel  to  the  azimuth  direction. 

Wave  number _ smooth ing  refers  to  ave’aging  neighboring  energy 

densities  of  the  transform  of  an  azimuthal  cut.  The  average  is 
over  energy  densities  at  equal  number  of  points  on  each  side  of 
the  given  wave  number. 

Another  wave-number  domain  averaging  is  averaging  of 
spectra  of  different  cuts.  In  this  method,  the  corresponding 


energy  densities  of  the  different  transforms  are  averaged.  This 
is  called  Wave  number  averaging. 


2.5.  Energy  Considerations 

The  Total  Energy  of  the  Signal  is  defined  as  E  | f ( n ) | 2 , 
where  f(n)  are  the  space  domain  sampled  signal.  The  Total  Energy 
of  the  Spectrum  is  taken  to  be  ^  x  | F ( k )  |  2  where  F(k)  is  the 
discrete  fourier  transform  of  f(n)  and  k1  includes  all  spectral 
points  except  k=0.  From  Parseval's  relation  for  discrete  transforms 

Y.  |  f  (  n  )  |  2  =  N  £  |  F  (  k  )  |  2  =  N(£!|F(k)|2  +  |  F  (  0  )  |  2  . 

The  total  signal  energy  remains  the  same  for  all  the  focuses  for 
a  given  range.  This  is  true  for  all  the  range  cuts  examined. 

2.6.  Opimalitv  of  Processing  In  Terms  of  Signal  To  Noise  Ratio 

Different  data  processing  methods  were  carried  out  to  reach 
a  method  that  generates  a  spectrum  w  i  t  li  a  high  signal  to  noise 
ratio-SNR.  The  different  processing  methods  tested  are  briefly 
d i s  c  u  s s  e  d  below. 

In  processing  trials  with  3  range  cuts,  the  cuts  are  taken 
at  ranges  =  64  ,  128  ,  192  (corresponding  t  n  incidence  angles  =  32.  5n, 

34.9,  3 7 . I r ,  respectively).  The  discrete  Fourier  Transforms, 

l)FT,  of  each  cut  is  taken.  The  spectral  energy  is  calculated 
and  smoothed  over  5  wave  numbers.  The  smoothed  spectra  are  then 
no mall  zed  and  averaged.  The  averaged  spectrum  is  plotted. 

For  processing  trials  with  5  range  cuts,  the  cuts  are  taken 
at  ranges  =  64,  96,  128,  160,  192  (corresponding  to  incidence 

angle  =  3  2.5  ,  3  3  .  7 r  ,  34.9’  36.0°,  37.1°,  respectively).  The 

procedure  mentioned  above  for  3  nit  s  is  then  repeated  for  the 
5  cuts. 
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For  processing  trials  with  5  range  cuts  and  no  area 
normalization,  five  cuts  are  taken  at  ranges  =  64,  96,  128,  160, 

192.  The  DFT's  of  each  cut  are  calculated.  The  energy  spectrum 
is  obtained  and  smoothed  over  5  wave  numbers.  The  smoothed 
spectra  of  the  5  cuts  are  then  averaged. 

In  processing  with  lateral  pixel  averages  (referred  to  as 
IA),  where  1A  =  40  at  range  128  implies  a  range  cut  at  128  that  is 
averaged  over  81  pixels  (128  -  40  pixels).  The  DFT  is  taken  for 
this  space  averaged  signal  and  the  energy  spectrum  then  calculated. 
Similarly,  for  IA  =  19,  the  procedure  implies  that  the  averaging 
is  done  over  39  pixels  (128  -  19).  For  I A  =  13,  the  same  procedure 

implies  averaging  over  31  pixels  (128  -  15). 

The  above  processing  methods  are  compared  in  terms  of 
SNR  values,  where  the  SNR  is  defined  as  the  ratio  of  the  energy 
density  at  the  dominant  wave  number  to  the  background  energy 
density.  A  sample  energy  spectrum  with  the  background  level  shown 
is  given  in  Figure  6.  The  SNR  values  for  different  processing 
trials  are  given  in  Table  5.  From  Table  5  it  is  seen  that  IA  =  15 

processing  trial  produces  the  largest  SNR  values  for  all  focus 

settings.  Hence,  this  averaging  trial  is  adopted  as  the  optimum 
processing  technique  for  all  the  SAR  data. 


0.0 


Figure  6. 


A  Spectrum  With  Noise  (background)  Level  Indicated 
Approximately  (SNRs7.1). 
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Table  5.  SNR  as  a  Function  of  Focus 
For  the  Different  Processing  Methods. 
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3 .  Results 

In  this  sot- Lion  tin*  results  from  tile  adopted  process  inn 
method  are  discussed.  The  dependence  o f  the  signal  energy  on 
incidence  angle  (range)  and  on  focusing  is  presented.  The 
dependence  o  f  the  dominant  wave  energy  density  on  focus  is  shown. 
The  de  pencil*  nee  of  Lite  energy  in  the  non -zero  part  of  the  wave 
number  spectrum  is  calculated  vs.  range  and  foe  us  selling. 


3.1.  Dependence  of  Energy  Density  at  Dominant  Wave  On  Focusing 

The  plot  of  the  energy  density  as  a  function  of  I 
is  given  in  Figure  7.  Here,  the  dominant  wave  number 
wave  that  is  140m  long.  Figure  7  shows  energy  density 
with  no  wave  number  smoothing  (N  =0),  and  with  wave  mi 
over  3  contiguous  wave  numbers  (N  =  3).  The  values  at 
I o  c  u s  corr e  spending  to  N  =  0  a  nd  N  =  3  ure  19.3  a  n d  I  0 

W h  e  n  defocused,  the  spectral  d e n  s  i  t  y  v a  1  u e  s  are  b  .  2  a  n 

respectively  corresponding  to  N  =  0  and  N  =  3.  The  local  ion  of 
the  peak  values  are  in  the  expected  region  according  to  the 
theories  presented  bv  .lain  (1979)  and  possible  llargi*  r  (1983). 
Symmetric  extension  of  the  curves  gives  a  central  peak  between 

fores  +3  and  +4  for  both  N  =  0  and  N  =3  curves,  shown  in  Figure 

7 .  These  results  indie  a  t  e  that  t  h e  opt  i mum  focus  setting  is 
between  .75c*  (focus  +3)  and  C  (focus  +4),  where  ('.  is  the  phase 
speed  of  the  dominant  ocean  wave. 
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The  I'eak-t  o-li,ise-Rat  io  (I*  HR)  Curves  are 
and  N  =  3  in  Figure  8.  The  I’ HR  is  defined  as 
energy  density  at  the  dominant  wave  number  to 
of  the  dominant  wave  number  for  the  do  focused 
corresponding  P  lilt's  are  3.2  1  and  2.14  for  N  = 
respectively.  The*  optimum  focus  is  approx  i  mat 
velocity.  For  comparison,  the  l’li  K  values  by  .1 
(1983)  are  shown  in  Figure  8. 
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Figure  7.  Energy  Density  at  X=140m  vs.  focus; 


6  t  =  34.9. 
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Figure  8.  PBR  at  A  *  140m  and  MARSEN  PBR  vs.  Focus.  (IA-15, 

Range= 128). 
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3.2.  Total  Energy  Dependence  on  I nc idence  Angle  and  Focus  Setting 

The  total  energy  decreases  mono t on ic a  1 1 y  with  increasing 
range  (increasing  incidence  angle-G^),  as  shown  in  Table  6. 

This  total  energy  is  insensitive  to  focusing.  If  the  D.C.  is 
removed  the  remaining  energy  at  each  focus  setting  exhibits  a 
local  maximum  at  range  =  128  (0^  =  34.9°),  as  shown  in  Figure  9. 
This  peak  is  sharpest  for  the  optimum  focus  setting  (+4). 

Without  D.C.,  the  remaining  energy  also  shows  a  dependence 
on  focus  setting.  This  dependence  is  plotted  for  the  cut  at 
range  =  128  (0^  =  34.9)  in  Figure  10.  This  Figure  shows  the  total 
energy  to  vary  smoothly  as  a  function  of  focus  setting,  and  to 
reach  a  peak  at  focus  setting  (+3).  The  equivalent  PBR  graph  is 
shown  in  Figure  11,  where  the  maximum  PBR  value  is  1.96.  Here, 

PBR  is  defined  as  the  ratio  of  the  total  energy  of  the  spectrum 
without  D.C.  to  the  total  energy  of  the  defocused  spectrum  (with¬ 
out  D .  C  .  )  . 
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Table  6.  Dependence  of  Total  Energy  (including 
D.C.)  on  Incidence  Angle. 
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Figure  9.  Total  Spectrum  Energy  vs.  Range  (14=15). 
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Figure  10.  Total  Spectral  Energy  vs.  Focus 
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.  Summary  and  Conclusions 


Different  processing  methods  have  been  tested  on  one  image 
data.  Processing  with  IA  =  15  (See  section  2.6.)  produced  the 
largest  signal  to  noise  ratio.  The  dominant  wave  length  is  found 
to  be  140  meters.  The  plot  of  the  energy  density  at  the  dominant 
wave  number  as  a  function  of  focus  setting  indicates  a  peak  energy 
density  value  at  a  speed  that  is  of  order  of  wave  phase  speed. 
Also,  the  non-D.C.  signal  energy  variation  with  focus  peaks  at 
focus  setting  of  the  order  of  wave  phase  speed. 
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APPENDIX  VII 


SAR  IMAGERY  SIMULATED  FROM  TWO-SCALE  RADAR  WAVE  PROBE  RETUREN 

William  J.  Plant 

Naval  Research  Laboratory 
Washington,  D.C.  20375 


The  coherent  ouput  from  NRL's  L-Band,  pulsed  radar  which 
was  mounted  on  the  NOSC  tower  during  the  TOWARD  experiment  was 
used  to  simulate  SAR  Imagery  of  the  ocean.  This  analog  output  of 
the  radar  was  beat  down  to  a  60  Hz  offset  frequency  and  recorded 
on  analog  FM  magnetic  tape  during  the  experiment.  In  the 
laboratory  it  was  digitized  at  a  rate  of  1024/T  Hz  where  T  is  the 
maximum  SAR  integration  time  for  which  simulation  is  desired. 
Digitized  data  were  stored  on  nine  track  magnetic  tape  and 
subsequently  read  into  the  simulation  program.  There  the 
simulation  proceeded  as  detailed  in  Plant  and  Keller  (JGR,  88, 
p.  9976,  1983).  Basically,  the  raw  data  were  both  averaged  after 

squaring  and  Fourier  transformed  over  an  integration  time  which 
could  be  T,  T/2,  or  T/4.  The  transform  was  performed  after 
clipping  the  data  and  multiplying  by  a  chirped,  offset  cosine 
signal  whose  chirp  frequency  was  related  to  the  focus  adjustment. 
The  power  spectrum  derived  from  this  transformed  signal  was  then 
convolved  with  the  squared  and  averaged  signal  which  had  been 
normalized  to  yield  cross  section.  This  process  gave  a  simulated 
value  for  the  spectrum  every  T/4  seconds  according  to  the 
relation 


1(0 


j  <XOO 


where  ^  *  was  determined  from  the  time  hy  multiplying  by  the  phase 
speed  of  the  dominant  ocean  wave  and  V  (  Xo  -  X  ) 
was  obtained  from  the  above-mentioned  power  spectrum 
Vt<r)  using  the  relation 


f  = 


V  ( x „  -  x ) 


7 r  /?< 


Note  that  only  simulations  of  az  muthal  cuts  through  a  SAR  image 
are  produced  by  the  above  procedure.  Subsequent  processing 
involved  obtaining  power  spectra  of  the  simulated  SAR  intensity 
images  and  averaging  several  of  these  spectra  together  to  form 
the  final  clutter-averaged  SAR  image.  Examples  of  the  results  of 
these  simulations  are  shown  in  the  following  figures  along  with 
power  spectra  of  the  cross  sections  obtained  during  the  process. 
All  simulations  shown  correspond  to  the  38,000  ft  flights  of  the 
SAR  on  October  31,  1984. 


Effect  of  Focus  on  Simulated  Image  Spectrum 
V-k  =  cos(ll°).  Leg  No.  1,  10/31/84,  Single 
Look,  N  =  27,  V/R  =  0.019  Hz.  T  =  2.2  see. 
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Figure  2  .  Effect  of  Focus  on  Simulated  Image  - 

V-it  =  cos (-  1  76  0  )  ,  Leg  No.  2,  10/31/84,  Singl 

Look,  N  =  46,  V/R  =  0.015  Hz,  T  =  2.6  sec. 
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Effect  of  Focus  on  Simulated  Image- 

V-k  =  cos(  158°),  Leg  No.  4,  10/31/84,  Single 

Look,  N  =  29,  V/R  =  0.015  Hz,  T  =  2.6  sec. 
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Figure  4.  Effect  of  Focus  on  Simulated  Image 

V-it  -  cos(-142°)  i  Leg  No.  7,  10/31/84,  Single- 

Look,  N  -  36,  V/R  =  0.015  Hz,  T  *  2.6  sec. 
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Figure  5 


Spectra  of  Cross-Section  at  Zero  Focus  with 
Conditions  Corresponding  to:  (a)  Figure  1, 
(b)  Figure  2,  (c)  Figure  3  and  (d)  Figure  4. 
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SAMPLE  PREDICTIONS  AND  SIMULATIONS  OF 
SAR  OCEAN  IMAGERY 


Robert  0.  Harger 

Electrical  Engineering  Department 
University  of  Maryland 
College  Park,  Maryland 


SUMMARY 


In  [1]  and  [2]  models  for  the  SAR  image  of  a  sea  surface  com 
posed  of  a  random  ensemble  of  small  gravity  waves  propagating  on, 
and  modified  by,  a  sinusoidal  long  wave  were  developed  from  funda 
mental  models  of  hydrodynamics,  electromagnetic  scattering,  and 
SAR  theory.  Generally  speaking,  the  hydrodynamic  and  electromag¬ 
netic  models  assume  "two-scales"  and  are  therefore  especially 
applicable  to  L-Band  systems  and  appropriate  seas-  as,  e.g.,  en¬ 
countered  during  the  TOWARD  experiment. 

The  detailed  analysis  of  [2]  yielded  a  fairly  complete  under 
standing  of  the  nature  of  these  models.  The  behavior  exhibits  a 
diversity:  a  recapitulation  will  not  be  attempted  here.  However, 
briefly,  it  was  shown  that  a  major  influence  on  the  SAR  image 
nature  was  the  "wavenumber  band  limiting"  effected  by  the  SAR 
system  on  the  phase  modulation  factor  due  to  the  round-trip  ray 
patch  to  the  large-scale  surface.  In  electrical  engineering  par¬ 
lance,  there  will  be  an  "output"  when  this  phase's  derivative-  a 
"local  wavenumber"-  lies  in  the  wavenumber  band  passed  by  the  SAR 
When  the  rate  of  change  is  sufficiently  large,  the  "output"-  if 
otherwise  unperturbed-  would  approximate  the  SAR  impulse  response 
The  "perturbations"  can  be  described  analytically  and  lead  to 
predictions  on  the  image  nature  and  details  such  as  the  "focus¬ 
ing  problem". 


A  principal  prediction  of  [2]  was  that,  under  appropriate 
conditions  on  SAR  and  long  wave  parameters,  "well-focussed"  ex¬ 
trema  can  occur  in  the  image  intensity  at  locations  corresponding 
to  the  crests  of  the  long  wave,  while  "dispersed  and  defocussed" 
contributions  occur  around  locations  corresponding  to  the  troughs 
of  the  long  wave.  In  particular,  the  average  (over  the  small 
scale  ensemble)  can  have  this  feature. 


The  following  is  an  excerpt  form  information  circulated  infor 
mally  within  the  TOWARD  group  in  summer,  1985,  and  contains  predic 
tions  of  certain  properties  of  the  SAR  images.  The  discussion 
assumes  that  the  long  wave  can  be  approximated,  at  least  locally, 
by  an  appropriately  chosen  sinusoid. 
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CHOICE  OF  RUNS  FOR  INITIAL  ANALYSIS 


The  runs  (2,6,11)  of  October  31,  1984,  were  chosen  for  an  Initial  "quick  look" 
because  (i)  there  is  evident  "wave-like"  imagery  of  good  quality,  especially  at  the 
38,000  ft.  altitude,  and  (il)  the  long  wave  is  azimuth-directed*  this  has  been  the 
more  "controversial"  situation  in  SAR  theory  and  practice  and  is  (there fore) the  most 
completely  elaborated  case  of  my  theory. 


Parameter  choices  for  analysis  and  simulations  -  The  following  choice  of  a  para¬ 
meter  set  was  made  for  this  preliminary  study.  ““ 

aircraft  parameters*  height*  38,000  ft.  =  14,140  m 

20,000  ft.  =  7,442  m 
7,000  ft.  =  2,605  m 
velocity*  231.66  in/sec 
"crab  angle"*  0 

longwave  parameters*  wavelength*  94  m 

amplitude*  0.435  m 

dispersion  relation*  gravity  wave,  18  m  depth 
direction*  "along  leg  1" 
shape*  sinusoidal 

small  scale  wave  parameters*  (Bragg  wavelength  =  28.83  m_1) 
nominal  phase  velocity*  0.5831  m/sec 

(wind  direction*  "along  leg  1") 

"white"  spectral  density  over  influential  wavenumbers 
"single  family"  with  "-a"  gravity  dispersion  relation 
no  current  effects 

SAR*  resolution  ("4- look")*  10  m  „  M 

spectral  weighting*  Gaussian  shaped,  truncated  at  3o  =  3/10  wave- 

number  bandwidth 
incidence  angle*  35 
beam  direction*  normal  to  flight  path 

By  way  of  further  comment,  the  longwave  wavelength  was  chosen  on  the  basis  of  the 
nominal  data  supplied,  the  approximate  spectrum  supplied,  and  checked  against  a  com¬ 
parison  of  SAR  imagery  and  bathymetry.  The  longwave  amplitude  was  chosen  so  that  its 
sinusoidal  form  would  have  the  same  rm’height  as  the  (average)  observed  long  wave* 
the  "significant  wave  height"  (nominal)  was  supplied  by  Scripps.  The  assumption  of 
the  lack  of  significant  current  was  based  on  current  information  supplied  by  Scripps. 
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NOMINAL  PREDICTIONS  OF  THEORY 

As  discussed  earlier  —  see,  e.g.,  my  hCMORANDUM  to  the  TOWARD  Data  Analysis 
Committee,  dated  June  5,  1985  —  when  the  longwave  amplitude  (A)  is  sufficiently 
large,  a  stationary  phase  approximation  usefully  describes  the  SAR  image  contri¬ 
butions  of  each  wavelet.  (The  required  condition  is  /KB>zlt  here  /KB  =  4.232|  on  the 
basis  of  experience  with  the  simulation  program  —  which  does  not  use  this  approxima¬ 
tion  —  the  condition  is  adequately  met). 

Such  a  wavelet’s  image  contribution,  recall,  is  described  by  a  doublet  (pair)  of 
SAR  azimuth  impulse  responses,  periodically  repeated  at  the  longwave  interval,  and 
translated  in  azimuth  by 


<  ]  k 
m*l0  ■"* 


VKmn> 

A 

K 


(1) 


and  degraded  by  the  quadratic  phase  "error" 


(2) 


6  *  ♦  1  for  one  member  of  the  doublet  and  6  =  -  1  for  the  other. 
Here,  with  little  error, 


oy(2k0sin«0) 

- 3T - 
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3.760  m,  Rn  =  2,605  m, 
10.743  m,  Rj;  =  7,442  m, 
20.412  m,  Rq  =  14,140  m, 


The  Influential  wavenumbers  of  the  small-scale  are  so  that 


I  6  \  ♦  (1  ♦  AB)  K2Q  «  1.314  m 
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Then,  in  (1),  for  the  three  altitudes  in  order  as  in  (3), 
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Therefore,  we  have  the  bound 
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6  =  -  1. 


Because  k  may  be  of  either  sign,  the  M spread"  of  the  (degraded)  SAR  impulse  loca¬ 
tions  os  ffwice  this  number.  AS  the  longwave  wavelength  is  94  m  and  the  SAR  resolu¬ 
tion  10  m,  we  see  thati 

1  -  When  6  =  -  1,  the  locations  are  fairly  well  clustered,  especially  so  at  the 

highest  altitude  and  less  well  at  the  lowest  altitude i 

2  -  When  6  =  +  1,  the  locations  are  rather  widely  spread  out,  especially  so  at 

the  highest  altitude. 

We  now  turn  to  an  analysis  of  the  effect  of  the  quadratic  phase  "error”.  This 
can  be  seen  by  calculating  the  resultant  rms  "resolution"  of  the  SAR  impulse* 
assuming  an  un-truncated  Gaussian  shape  for  the  azimuth  transfer  function,  we  get  the 
approximation 
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Here,  then, 
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12 . 28m 
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1 0 . 0@m 
10.17m 
.1 1 . 02m 


,  «  =  -  1 

We  see  that  both  members  of  the  doublet  are  degraded  rather  little  at  all  the  altitu¬ 
des,  for  this  particular  data  set. 

We  can  now  make  some  predictions  of  the  image  nature  based  on  the  above  simple 
analysis. 


The  image  nature  is  primarily  affected  by  the  translations 
( wavenumber-dependent)  of  the  individual  imaae  contributor 


(wavenumber-dependent)  of  the  individual  image  contributions  of  each  wave- 
letj  it  is  affected  rather  little  by  the  degradations  of  the  doublets. 


The  image  will  have  local  maxima,  spaced  at  the  longwave  wavelength. 
Accounting  for  the  relatively  small  translations  of  these  impulse  respon¬ 
ses,  their  "resolution",  or  rms  spread,  will  be  estimated  as  about 


36  m,  Rn  =  7,000  ft.  =  2, 
'  24  m,  R~  s  20,000  ft.  =  7, 
13  m,  Rq  =  35,000  ft.  =  14, 


605  m, 
442  m, 
140  m, 


-  The  "background"  against  which  these  local  maxima  appear  will  be  of 
increasing  spread  —  hence  (to  a  limit)  of  reduced  level,  as  the  height 
increases. 

-  Therefore,  the  "contrast  ratio"  should  increase  with  increasing  height. 


Thus,  we  expect  to  see  at  all  three  heights,  local  (maxima,  spaced  at  the 
lowes tVheightlen9th *  WUh  contrast  ratio  best  at  the  latest  height,  worst  at  the 


The  optimum  value  of  the  focus  correction  parameter  C 
MEMORANDUM  —  is,  here, 


—  see  the  above-cited 


24.1  *  10‘ 


=  -10.4*  10'*  +  8.44  *  10-2  ,  R~  = 

4.44  *  10-2  .  R_  = 


2,605  m, 
7,442  m, 
14,140  m, 
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In  the  TOWARD  experiment  the  observed  long  wave  spectral 
density  was  apparently  not  especially  narrow  in  its  support, 
giving  impetus  to  extending  the  theory  of  [2]  in  this  regard: 
this  has  now  been  done  under  appropriate  approximations.  A  basic 
difficulty  is  that  an  adequate  theory  for  the  modification  of 
short  waves  on  a  suitably  general  long  wave  is  even  now  under 
development.  (See  Longue t -H igg ins ,  Appendix  IV  herein). 

If  the  bandwidth  of  the  long  wave  ensemble  is  not  unduly 
broad,  its  sample  functions  will  have  a  quasi-sinusoidal  nature 
and  we  expect  that,  locally,  the  main  features  of  the  above 
analysis  carry  over:  this  is  borne  out  by  the  simulations. 

Some  examples  of  simulation  results  will  be  discussed  below 
for  a  set  of  parameters  that  include  an  assumed  5m  resolution. 

The  JPL  SAR  employed-  at  least  some  of  the  time- 
a  very  widebeamed  antenna  having  an  intrinsic  resolution  about 
ten  times  finer,  while  the  imagery  distributed  to  the  investi¬ 
gators  was  degraded  to  various  coarser  resolutions-  e.g.,  11  m. 

The  effects  of  motion  in  such  data  can  be  approximately  equiva¬ 
lent  to  those  of  a  SAR  with  such  a  resolution-  e.g.,  11m. 

A  more  precise  simulation  would  likely  require  the  program  to 
be  transported  to  a  computer  of  significantly  greater  capacity. 

(A  VAX  785  is  currently  being  used). 
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The  pricipal  parameters  assumed  in  the  simulations  cited 
below  are  listed  here. 


SAR  resolution:  5  m 

SAR  slant  range:  14,000  m 

SAR  RF  wavelength:  0.25  m  (L-band) 

SAR  incidence  angle:  35  deg 

Long  wave:  (a)  sinusoidal,  wavelength  100  m, 

amp litude  0.5  m 

(b)  random  with  wavenumber  spectral  density 
of  Scripps  for  Oct.  31,  1984 
water  depth:  18  m 

Small  scale  ensemble:  random  with  "white"  spectral 
density,  one  "sigma  family" 

Simulation  sample  distance:  1  m 
Number  of  object  samples:  1024 
Number  of  image  samples:  338 


In  Figure  la  is  shown  the  simulated  image  for  the  sinusoidal 
long  wave  and  a  small-scale  wavelet  consisting  of  a  single  sinusoid 
at  the  "Bragg-matched"  wavenumber.  Its  nature  is  fully  explained 
by  the  analysis  of  [2]-  recall  the  image  extrema  correspond  to  the 
crests  of  the  long  wave.  The  corresponding  spectrum  (Fourier  trans¬ 
form)  is  shown  in  Figure  lb:  of  course  it  has  strong  "lines"  as  a 
result  of  the  apparent  periodicities  in  the  image.  (Note  that  the 
harmonic  relation  to  the  long  wave  becomes  more  pronounced  at  the 
higher  wavenumbers.  Also  recall  that  a  long  wave  traveling  with  the 
SAR  has  wavelength  its  increased-  here  a  few  percent). 
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In  Figure  2a  is  shown  the  sample  average  of  about  two  dozen 
images  for  the  same  sinusoidal  long  wave  but  with  a  full  spectrum 
of  small-scale  wavelets.  Again,  the  analytic  predictions  of  [2] 
are  supported  by  the  numerical  simulation.  The  corresponding 
spectrum  is  shown  in  Figure  2b:  the  "lines"  due  to  the  image  period¬ 
icities  are  still  clear,  though  not  as  abundant. 

In  Figure  3a  is  shown  a  segement  of  a  sample  function,  at  the 
mean  observation  instant,  of  the  long  wave  random  ensemble.  This 
(Scripps)  spectrum  has  considerable  support  and,  correspondingly, 
the  sample  functions  fluctuate  more  rapidly  than  a  single  sinusoid 
at  the  dominant  wavenumber.  Based  on  the  analysis  of  [2],  we  expect 
image  maxima,  generally,  corresponding  to  the  crests  of  the  long  wave 
sample  function. 

In  Figure  3b  this  expectation  is  borne  out.  The  theory  of  [2] 
predicts  an  image  translation  of,  here,  about  21  m  due  to  the  mean 
effect  of  the  small-scale  wave  motion:  if  3b  is  so  translated,  its 
extrema  lie  rather  as  predicted. 

The  corresponding  spectrum  is  shown  in  Figure  3c:  the  "lines" 
appearing  are,  at  the  higher  wavenumbers,  generally,  harmonically 
related  to  a  wavelength  about  100  m-  approximately  the  dominant 
wavenumber.  However,  at  the  lower  wavenumbers  there  is  strong 
spectral  content  attributable  to  the  longwave  and  image  fluctuation 
which  is  a  fraction  of  the  long  wave  period. 

A  (limited  as  yet)  number  of  similar  simulations  indicate  that, 
rather  as  would  be  expected,  the  image  nature  will  strongly  depend 
on  the  particular  long  wave  sample  function  "drawn  by  nature"  from 
the  ensemble  charaterized  by  the  spectral  density  measured  by 
Scripps.  Such  a  situation  calls  for  a  suitable  averaging  process 
over  the  long  wave  ensemble. 
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In  Figure  4  is  shown  a  "standard"  ([3],  p.  556)  spectral  density 
estimate  using  24  complex,  simulated  SAR  images.  on  each  run  both  the 
long  wave  and  the  small  scale  were  drawn  randomly. 

We  may  compare  this  spectral  density  estimate  based  on  the  simu¬ 
lation  with  the  actual  SAR  image  data.  Presently  at  hand  is  a  spec¬ 
tral  density  estimate  obtained  by  the  same  "standard"  method-  but 
using  100  adjacent  range  cuts:  therefore,  because  of  the  nearly 
azimuth-directed  nature  of  the  long  wave,  the  long  wave  samples  are 
rather  correlated.  Such  a  spectral  density  estimate  is  shown  in 
Figure  5. 

Comparing  Figures  4  and  5  we  see  that  the  SAR  image  spectrum 
is  "full"  to  wavenumbers  cutoff  by  the  bandwidth  of  the  SAR  system. 
Also,  comparing  to  the  long  wave  spectral  density  measured  by  Scripps^ 
appearing  elsewhere  in  this  report,  we  see  that  it  does  not  resemble 
at  all  either  the  "actual"  or  "simulated"  SAR  spectral  densities. 

While  subject  to  refinement,  we  expect  the  above  conclusion  to  hold 
and  is  a  very  important  observation  since  many  have  mistakenly 
believed  that  they  are  similar. 

As  a  final  simulation  example,  a  SAR  image  of  344  m  (azimuth) 
by  210  m  (range)  is  shown  in  Figure  6.  (The  aspect  ratio  is  incor¬ 
rect,  being  taken  "as  is"  from  the  printer).  The  (randomly  occur¬ 
ring)  maxima  attributable  to  the  long  wave  are  evident. 

This  simulation  capability  could  be  a  powerful  assist  in  under¬ 
standing  SAR  imagery  of  ocean  phenomenon  in  various  situations  and 
applications . 
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